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ALGORITMO F.F.T. PARA DETERMINAÇÃO 
DA TRANSFORMADA DE FOURIER -— SUA 
PROGRAMAÇÃO EM FORTRAN 


RESUMO 


Além de uma explicação passo por passo da progra- 
mação, em Fortran, do algoritmo de F.F T. também, 
se apresentam alguns resultados obtidos com o programa 
proposto que permitem julgar da excelência do método, 


1. INTRODUÇÃO 


Desde há muito que a transformada de Fourier 
é utilizada para caracterizar os sistemas lineares 
e para identificar as componentes das várias 
frequências existentes numa forma de onda 
continua. Ao recorrer-se a um computador 
digital a análise tem que se resumir a um 
conjunto finito de pontos. É necessário proceder 
a uma amostragem da forma de onda ou à aná- 
lise por pontos do sistema linear recorrendo 
então à versão discreta da transformada de 
Fourier, que designaremos pelas iniciais D.F.T. 
(«Discrete Fourier Transform») tal como na 
literatura de língua inglêsa. 

O cálculo da transformada discreta de Fourier 
pode fazer-se através de um algoritmo, conhecido 
pelas iniciais F.F.T. («Fast Fourier Transform»), 
cuja eficiência se encontra amplamente discutida 
na literatura da especialidade, quer em relação 
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SYNOPSIS 


It presents a thoroughly explanation of the Fortran 
trogramming for the F.F.T. algorithm. In order to 
evaluate its performance it still presents some results 
obtained by the proposed Fortran program, 


ao tempo de cálculo, quer em relação ao erro 
proveniente dos arredondamentos. 

O interesse da transformada de Fourier e, 
portanto, do algoritmo F.F.T. para a sua deter- 
minação é enorme na Análise de Sinais, assunto 
a que se dedica um dos autores no Centro 
de Estudos de Electrónica da Comissão de 
Estudos de Energia Nuclear (Instituto de Alta 
Cultura). Das aplicações destacam-se a determi- 
nação do espectro de potência, a determinação 
de integrais de Fourier, séries de Fourier e inte- 
grais de convolução, a simulação de filtros e a codi- 
ficação da voz. 

Foram estas as razões que levaram à escolha 
deste assunto para um trabalho que, ao abrigo 
da rubrica de Actividades do III Plano de Fomento, 
se destina a estudantes que desejem aperfeiçoar 
e desenvolver os conhecimentos ministrados na 
cadeira de «Elementos de Análise Numérica» do 
Instituto Superior Técnico. Para um melhor apro- 


veitamento recorreu-se a estudantes que haviam 
cursado a cadeira, acima mencionada, no ano 
lectivo anterior. A utilização de computadores, 
em geral, e de programação Fortran, em particular, 
foram os assuntos dados na cadeira e que os 
estudantes visaram ao colaborar no trabalho e 
neste artigo. 

A programação do algoritimo F.F. T. foi rea- 
lizada tendo em mente o computador IBM 1620 
cuja possibilidade de utilização agradecemos à 
Companhia IBM Portuguesa e dentro da com- 
panhia ao serviço de cálculo científico, pois 
desta utilização depende em grande parte a rea- 
lização do objectivo deste trabalho. O programa 
Fortran obtido, ainda que tendo em vista o tipo 
e a dimensão do computador à nossa disposição 
é imediatamente adaptável a computadores de 
maior capacidade de memória, sujeitos apenas à 
restrição de utilizarem o Fortran. Além disso, 
mesmo para o IBM 1620, através de uma sub- 
divisão em vários programas (sequencialmente 
submetidos ao computador) foi possível estender 
o algoritmo para um número de pontos superior. 
Antevê-se mesmo a possibilidade de uma con- 
densação ainda maior do programa que só não 
se realizou por carência de tempo, uma vez que 
este trabalho foi realizado ao abrigo do III Plano 
de Fomento para 1968 e, portanto, de duração 
limitada. 


2. PROGRAMAÇÃO 


Na programação de um método de solução de 
um problema é necessário antes de mais analisar 
a teoria na qual se baseia o algoritmo em questão. 
Por isso começaremos por apresentar a teoria do 
algoritmo F.F.T.. 


2.1 Teoria do algoritmo F.F.T. 


Vamos ver o que vem a ser a D.FT. de 
uma função discreta. Sendo T o intervalo em 
que está definida a função discreta 


X (tr)=X, com ks:(0,1,2,...,(N—-1)) 
cuja transformada discreta de Fourier é 
V(g)=V (N-1)) 


então verificam-se as seguintes expressões que 
estabelecem as relações entre os pares de Fourier, 
no caso discreto, 


com selo, - E, TD 


N—1 


Vy=4t 3) X exp(—j2z7ks/N) 
k=o 
(1) 
gm20,1,2, «os tN=T) 
N-1 
Xe, = df >, V. exp(+4 j27skN) 
(2) 
k=0,1,2 :coAtN=T) 
em que 
f,=sáf 
t,=kA4t 
At = T/N 
Af= 1/N 


Vamos concentrar-nos no cálculo da transfor- 
mada discreta directa de Fourier, isto é, no cálculo 
de V, a partir de X,. Para maior comodidade 
utiliza-se At=-1 o que equivale a dizer que o 
intervalo T, nas unidades que vamos utilizar, 
mede exactamente N unidades. 

Analisemos o caso de N ser par; se fizermos 


| nd. 
Zm = Xom+1 (3) 
m=0,1,2, (5-1) 
2 / 
então (1) virá 
N-1 
Mg X. exp(—j2=ks/N) = 
k = 0 
N/2— | 
= [Ym exp(—já=ms N) + (4) 
m==o 
fo La expl—jzr(2m +1)sN]! 
8=0,1,2,...,(N—1) 
ou 
N2—1 
V;= 5 [Ye exp(—já 7 ks/N)) + 
k=o 
+ exp (— jz ms N) 
Ng-1 
> (Z, exp (— jázks/N)! (5) 
k=o0 
$==0,1,2,...,(N-—1) 
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Fazendo 


N2 1 
U,= bz Y. exp(—j4 -ks'N) 


k = o 


N2—1 (6) 
T,= > Z, exp(—j4ázksN) 


k=0 
em que s varia entre zero e (N—1), em prin- 
cípio. Dado que U,, N =U,eT,, N =T, 
2 2 
os únicos valores distintos são parast 10,1,2,..., 


N 
(5-1)| ; por esta razão podemos dividir a 
/ 


expressão (5) em duas 


V=U,+T, exp (— j2 7 s/N) 


Vad Ur exp | —j2= (s+ 2) NI] 
2 
$=0,1,2,..., (41) (7) 
2 


Finalmente para simplificar a notação pode 
escrever-se 


W(s) = exp(—j2=s/N) (8) 


pelo que se obtém 


V=U,4-T,W(s) 


Vis N=U—-TW(s) (9) 


2 
N 
se [0,1,2,... ( 2 1)] 
2 


s + 


Nisto 
N 
expl-iz=(s d E /N|] = exp (—izrs/N— 


A. lh EM e a 
—j22-0) exp(—-j2=s/N) exp(—j=)= 


/ 


= —exp(—j27=sN) (10) 
Ou então também podemos escrever 
V=U,+T, W(s) 
N 
Gs = UT W(s+5) (1) 
2 
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com se[0,1,2,..., (5-1) ) . 


Vamos ver um exemplo para N ==8 em que 
recorremos na Fig. 1 a uma exposição à base 
de grafos. Desta figura se conclui que transfor- 
mámos o cálculo de uma D.F.T. (N = 8) em duas 
D.F.T. (N==4) e num conjunto de operações que 
são expressas pelas equações (11). 


D.F.7 


[N =4) 


DF 1 
(N=4) 


Fig. 1 


Observando a expressão (6) conclui-se que 
existe uma analogia entre U, e T, e comparando 
as expressões (5) e (6) verifica-se uma seme- 
lhança daqueles com V.. Se existir um número 
par de valores de U, e de T, aquela semelhança 
vai-nos permitir simplificar o cálculo uma vez 
que U, e T, são transformadas discretas de 
Fourier (neste caso dos Y e dos Z). Quer dizer 
que podemos repetir para U, e para T, as opera- 
ções efectuadas anteriormente sobre V.. Nesse 
caso U, e T, obter-se-ão pelas expressões 


U,=R,+5S, W(2s) 


UN =R,+S,W] 2 (+ 5] 
4 /. 


s=P +Q, Wí(2s) 
| mA 
A N=P,+Q,W] 2 (s+4)| 


se/0, 1, Bee (a 1)| 


4 
Se for possível proceder sucessivamente de 
igual modo chegaremos ao cálculo de um conjunto 


de transformadas de Fourier para (N==2) para 
as quais os valores que entram nos cálculos são 
os próprios X, (valores da função dada). A possi- 
bilidade de subdividir sucessivamente por 2 os 
grupos de valores obtidos exige que N (número 
total de valores dados) seja potência inteira e 
positiva de 2. 


Outros algoritmos F. F. T. podem basear-se em 
potências de 3, 5 ou outro número primo qual- 
quer; todavia o cálculo simultâneamente mais 
simples e económico é aquele em que N é potên- 
cia de 2. O algoritmo F. F.T. é ainda aplicável 
quando N for um produto de vários factores 
primos mas então não existe uma repetição do 
cálculo, mas simplesmente uma repetição de 
conceitos ao longo do cálculo, pelo que este se 
torna muito mais complicado. Se N for um 
número primo o algoritmo F.F.T. e o cálculo 
directo coincidem. 

Em geral escolhe-se N potência de 2 pela sim- 
plicidade que traz ao algoritmo e à sua progra- 
mação e é essa também a nossa atitude. No caso 
de (N==8) o conjunto das transformações pode 
esquematizar-se no grafo da Fig. 2. 


a: 4 yo 
: A + " nx a 4 x Pay ) 
E A EA / 
(x) = 5.) Au) x =(v,) 
= “o E N Y w » 
(15) SAS) 4 » E X +v,) 
E A, X, 
| ) SF) riso AD, Sé x ATA v.) 
d x X N 
N N 
(45) WA) - reg ) Ko —Age( vo) 
- N N 
(M,) = 0.) >< - el 1 ) , x ra Va) 
(o) ma) a) a "3; 
Fig. 2 


Podemos fazer uma extensão das expressões 
(9) para o cálculo total de uma F.F.T. com 
N = 2*É tendo em consideração o grafo da Fig. 2 
e indicando pela variável i as colunas (a partir 
da direita). Resultará 


Vi-D=VMA+V (O) Wis 27?) 
Vesgli-D=V (DV, ()W(s27?) (13) 


sendo V,(i) o valor que ocupa a coluna i e a 
linha r no grafo correspondente ao da Fig. 2. 
Deve-se ter 


N 


a $i 


KK = log, N 
ic(2,3,4, 064 KKK 4 1)) 
re(1,2,3,...,9] 


=[—1l 
u=rT 
v=r+0 


Esta variação de r só é válida para uma parte 
dos N valores que é necessário calcular em cada 
coluna de valores. Assim na coluna de índice i 
necessitamos de 2'”* variações de r cuja forma 
geral será 


rej2(k—1)9+1,2(k— 1)9+2,...,(2k—1)5) 


com k:[1,2,3,...,27?). O índice r utilizado 
nestas expressões tem uma variação que começa 
com a unidade, o que é vantajoso pois também 
em Fortran os indices começam com a unidade, 
não podendo começar em zero tal como acontecia 
com o índice s. 

Por conveniência futura vamos indicar as 
colunas de valores a calcular pelo indice j que, 
contrariamente ao índice i, cresce da esquerda 
para a direita e definir uma nova variável 
A, ()=—V, (KK —j + 2). Então as expressões 
(13) podem reduzir-se a 


A(D)=A(-D+A,G-1) W(s 2847)) 


(1, 2,30 con KK) 


N 
k € da Di Meyers À 
Let t, 2, 3;.00, 2] (14) 


q L «1 
r=(k-D27+L 


Larg =>[10 =! 
v=[+2'" 
o uiaÃ u=r—2'"' 
Ls? ==> 
Vv=r 


As expressões acima são a base do algoritmo 
F.F.T. pelo que encerramos as considerações de 
ordem teórica. 


2.2 Execução da programação 


Uma vez indicada a estrutura teórica do algo- 
ritmo vamos delinear o programa. 

Pretende-se um programa que aceite dados 
em número variável com as restrições de ser 
potência de 2 e de ser inferior a um determinado 
número imposto pela dimensão da memória do 
computador a usar. 

No programa que nos propomos fazer, quer os 
dados, quer os resultados e, portanto, todos os 
valores intermédios do cálculo são números com- 
plexos. Possuem parte real e parte imaginária 
que em Fortran identificaremos por um índice: 
1 para a parte real e 2 para a parte imaginária. 
Todas as quantidades acima referidas serão, pois, 
variáveis indexadas (ou dimensionadas) em For- 
tran. 

Para se compreender as preocupações a que a 
seguir nos vamos referir será conveniente afirmar 
que as transformadas discretas de Fourier para 
serem de marcada utilidade devem poder utilizar 
um grande número de pontos. Será, pois, com- 
preensível que se procure um programa econó- 
mico em ocupação de memória do computador 
para que a restante parte da memória possa ser 
utilizada no sentido de aumentar o número de 
pontos iniciais do algoritmo F.F.T.. 

Ora, em princípio, com N = 2*% é necessário 
calcular sucessivamente KK + 1 conjuntos de N 
valores que já sabemos serem complexos e que 
ocupariam 2N (KK +11) posições de memória. 
Como, porém, os cálculos indicados pela expres- 
são (14) se efectuam entre um conjunto de valo- 
res e aquele que se lhe segue é possível reduzir 
os KK + 1 conjuntos de valores a 2 conjuntos 
de N valores complexos. Ocupam desta feita 4N 
posições de memória o que representa uma 
economia importante na ocupação da memória 
do computador. 

Iremos identificar estes dois conjuntos de 
valores, em Fortran, pelas variáveis indexadas A 
e X e atribui-se a A os valores do 1.º, do 3.º, 
do 5.º, etc., conjunto de valores e a X os valores 
do 2.º, do 4.º, do 6.º, etc. conjunto de valores. 
Logo calcula-se X a partir de A e depois A a 
partir de X e assim sucessivamente até atingir- 
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mos o conjunto de valores finais que será À ou 
X conforme KK seja par ou ímpar, respectiva- 
mente. 

Para identificar cada valor dentro do conjunto 
representado por A ou por X vamos utilizar um 
primeiro índice que irá dela N para além do 
segundo índice que identifica a parte real e a 
parte imaginária. Assim teremos 


A(LK) ou X(LK) 


com 


1:(1,2,3,..., N) 
Ke(1,2) 


O valor de W(s 2** —)y também é um valor 
complexo que representaremos, em Fortran, pela 
variável indexada W, cujo índice identifica se se 
trata da parte real ou da parte imaginária. 
Como a expressão de W(s 2*K—1) é dada 
por (8) a fórmula de Euler permite determinar 
como parte real e como parte imaginária, respec- 
tivamente, (*) 


W(1) > cos(— 2758 2*K—1/N) 
W(2) -— sin(— 275 6-1 nys 


Como vimos, a obtenção dos resultados pela 
ordem natural exige que os dados sejam orde- 
nados de uma forma especial. Esta ordenação 
será efectuada do conjunto X para o conjunto A, 
com economia de 2 N posições de memória, 
pois, não é necessário um novo conjunto para 
os dados iniciais. Desta forma o conjunto 
X é que recebe os dados iniciais a introduzir no 
computador como representando os pontos da 
função a transformar pelo algoritmo de F.F.T.. 

Como hipótese vamos estabelecer a forma 
como são introduzidos no computador o valor 
de N e os valores do conjunto X. N será intro- 
duzido através de um cartão perfurado e obede- 
cendo ao formato (14); para o conjunto X 
utilizaremos cartões perfurados que se encontram 
por ordem e em que existe um ponto em cada 


(*) Utilizaremos o símbolo <«--—» para indicar a equi- 
valência entre expressões ou variáveis, em lingua- 
gem Fortran, e as expressões ou variáveis, em 
linguagem usual 


cartão (1.º a parte real e depois a parte imagi- 
nária) segundo o formato (6 X,E 14.8,1X,E 14.8). 

Também estabeleceremos desde já o aspecto 
da saída dos resultados: por impressão numa 
folha de papel, em forma de tabela que indique 
um índice de ordem do argumento, a parte real 
e a parte imaginária da transformada. O formato 
a utilizar deve ser (5 X,15,10X,E 12. 6,10X, 
E 12.6) e a tabela deve ser antecedida pela 
captação 


J1 REAL IMAGINARIA 


segundo o formato( ///7 X,2HJ1,15X,4HREAL, 
15X, 10H IMAGINARIA /). 

Para indicar ao compilador Fortran que o pro- 
grama está terminado colocar-se-á em último 
lugar a instrução END. 

Dado que o computador IBM 1620 posto à 
nossa disposição para o trabalho tem uma capa- 
cidade limitada (20 K) o número de pontos máximo 
não pode exceder 128, uma vez que tem de ser 
potência de 2. Assim podemos escrever a instru- 
ção de especificação 


DIMENSION A (128,2), X (128,2), W (2) 


que indica a dimensão máxima das variáveis inde- 
xadas que entram no programa: A, X, W. 

Caso este programa se queira utilizar noutro 
computador “com maior capacidade de memória 
basta alterar a instrução DIMENSION indicando 
em vez do número 128 a potência de 2 que essa 
capacidade permitir. 

A elaboração do pormenor do programa vai 
ser apresentado separado em várias fases. 


2.2.1 


Nesta primeira fase analisar-se-á a introdução 
dos dados. Os dados a introduzir são o valor de 
N (número total de pontos dados) e os valores 
da parte real e da parte imaginária de cada um 
desses pontos. Como se utilizam cartões perfu- 
rados teremos em Fortran 


999 READ 11, N 
DO 121=1,N 
12 READ 13, X(1,1), X (1,2) 


em que se recorreu a um método iterativo para 
a leitura dos valores da parte real e da parte 
imaginária que existem em cada cartão. 


2.2.2 


Se examinarmos a Fig. 2 verificamos que, para 
obtermos a transformada de uma função dada 
por N pontos X,, X,, X., ..., Xy de modo a que 
os resultados venham segundo a ordem natural, 
é neçessário dar uma ordenação especial aos 
dados. Nesta segunda fase discutiremos a exe- 
cução do programa relativa à ordenação dos 
dados. 

Da análise da expressão (3) conclui-se que 
agrupámos os valores de índice par no conjunto 
dos U e os de índice impar no conjunto dos T. 
Se os colocarmos, tal como no grafo da Fig. 2, o 
conjunto dos U é seguido do conjunto dos T. 
No passo anterior é necessário operar em cada 
metade como anteriormente, isto é, separando os 
pares dos ímpares e colocando aqueles antes 
destes. 

Vamos agora supor que o índice de ordem dos 
valores V, se encontra escrito na forma binária, 
isto é, 


s=(Dex Du — 9 da , b,, bo) = 
= be 210 ba eos 
+b,2+b, (16) 


em que para vicío,1,2,...,(KK—1)Ise tem 
b,e[0,1] 


Repare-se que a divisão em índices pares e 
índices ímpares, se pode fazer atendendo que em 
escrita binária os primeiros têm b,=0 e os 
segundos têm b, 1. Ao fazer-se o mesmo em 
relação a cada um dos grupos assim obtidos 
separam-se os valores que têm b —O dos que 
têm b,=1, isto é, tendo em atenção o segundo 
símbolo a contar da direita. E dentro de cada 
grupo colocam-se em primeiro lugar aqueles em 
que b,==0. Depois, separar-se-iam os que têm 
b,=0O daqueles em que b,=1 e assim sucessi- 
vamente. 

Pode concluir-se que os valores iniciais se 
encontram na posição que corresponde a escrever 
o valor do índice em escrita binária e a inverter 
a ordem dos algarismos binários. 

Há, pelo menos, duas formas possíveis de pro- 
ceder a esta ordenação. 

A primeira corresponde a tomar o valor de 
ordem da posição a preencher, escrevê-lo em 


TECNICA Nº 391 


notação binária, inverter a ordem dos símbolos 
binários e determinar o valor que lhe corresponde; 
este valor indica o índice. Se tal como em For- 
tran os índices e as ordens começarem pela uni- 
dade e não por zero teríamos para N=8 o 
seguinte esquema: 


J=1;j=0j=(000)-->i=(000)i=01=0 
J=2j=1j=(001)=-> i=(100) i=41=5 
J=3j=2j=(010)-->i=(010)i=21=3 
J=4;j=3j=(011) => i=(110)i=61=7 
J=5 j=4j=(100)=>i=(001)i=11=2 
=6j=5j=(101)->i=(101)i=51=6 
J=7j=6j=(10)=->i=(0)i=31=4 
J=8j=7j=(11)=>i=(n)i=71=8 
(17) 


A segunda forma é através de um método de 
recorrência que determina o índice i, para a posi- 
ção de uma dada ordem j a partir do índice i,., 
conhecido para a posição de ordem imediatamente 
anterior j—1; será necessário determinar q tal 
que 

j=i, 1,19 (18) 

É este caminho que seguiremos, tanto mais que 
vo só depende, quando j se encontra em escrita 
binária, da posição do primeiro algarismo dife- 
rente de zero (a partir da direita). Pode verifi- 
car-se que o número de zeros que antecede este 
primeiro algarismo diferente de zero varia entre 
zero e KK — 2. 

Tal como já fizemos vamos designar por J a 
ordem da posição a preencher e por I o índice 
que lhe corresponde, tanto mais que estas desig- 
nações correspondem exactamente às que utili- 
zaremos no programa Fortran. 

Como se trata de um método de recorrência é 
necessário indicar as condições iniciais que neste 
caso se resumem a afirmar qual é o índice que 
corresponde à primeira posição J==1; do que 
foi dito se conclui que o índice é 1 


1 = 2 (19) 


Em Fortran resulta 


[| =1 
A (1,1) = X (1,1) 
A (1,2) =X (1,2) 
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A partir da posição de ordem 2 até à posição 
de ordem N vai repetir-se um conjunto de ope- 
rações que em Fortran sabemos poder indicar 
através de uma instrução 


DO 2] =2,N 


As operações de que falámos acima compreendem, 
em primeiro lugar, uma verificação referente 
à ordem J em que se procura saber qual a potên- 
cia de 2 de quej=]J—1 é um múltiplo ímpar. 
Desta forma identificamos a posição do primeiro 
algarismo diferente de zero, a contar da direita, 
existente na escrita binária de j. 

Com esse fim faz-se a variável K2 percorrer 
as potências de 2 com início em 2' 


KK =1 
1K2=2.**KK + ..01 


K2 A 
então e também percorre as potências de 2, 


mas agora com início em 2º. Se designarmos por 
L (tal como em Fortran) a diferença entre (] — 1) 


K2 
e(2KA-—1) Ei vamos verificar se L se torna 


igual a zero quando 2 KA — 1 percorre os núme- 
ros ímpares. Enquanto L for positivo é necessário 
passar ao número ímpar seguinte, mas se L se 
tornar negativo sem passar por zero, então J—1 
não é múltiplo ímpar dessa potência de 2 e 
passa-se à potência seguinte recomeçando com 
KA = 1. 


Mi=N/K2+1 
DO 3 KA = MI 
L=]-K2*KA+K2/2—1 
IF (L) 5,4,3 

3 CONTINUE 

5 KK=KK +1 
GO T0 1 


Se verificarmos que ] — 1 é múltiplo ímpar da 


potência de 2 dada por LÃ então sabe-se que 


após a inversão da notação binária se tem um 
conjunto de algarismos binários zero à esquerda, 
em que o primeiro algarismo diferente de zero 


a partir da esquerda é dado por — a dividir por 


K 2 
—— , ou seja —. 

2 K2 

Consideremos, agora, a ordem da posição prece- 
dente j—-1; quando em escrita binária, possui, 
a partir da direita, algarimos 1 onde j possui 
algarismos O e um algarismo zero na posição 
correspondente ao primeiro algarismo 1 a partir 
da direita em j. À esquerda deste algarismo ambas 
as escritas binárias coincidem. Após a inversão 
da notação binária para ambos os valores, a dife- 
rença será fácil de obter : 


N N 9 N 
== —- 2-2 — — 
K 2 K 2 K 2 
bp N (20) 
— 1. — 2 K2 
N N 
bH pc a 
em que 2 K2 e, portanto 2 K> a 
Abreviando pode-se escrever 
i—i, ,=a;—a, (21) 
em que 
N 
a, = 
K 2 
N N 2 
pego E duto 
K 2 K 2 2 


Como a, é a soma de uma progressão geomé- 
trica de razão 2, então 


N N 
2 xa ê a No (23) 
a, = emo: qe (NS e 
z L— 2 K 2 


Incluindo este resultado na expressão (21) 
obtém-se 


N N 
list sN+hã--= 
K 2 K 2 
N (24) 
=1,—(K2—3 
= o 


ou em Fortran, em que designamos por I qualquer 


dos i, 


41=1I-(K2 — 3) *(NK2) 


Determinado o índice Í que corresponde à 
posição ] resta indicar a arrumação dos valores 
iniciais X no conjunto À 


A (5,2) = X (1,2) 
2 A(,1)=X (11) 


Quando estas operações forem realizadas para 
todos os valores de ]J desde 2 até N tal como 
está indicado na instrução DO 2] == 2,N, então 
os valores iniciais encontram-se ordenados no 
conjunto A conforme convém para prosseguir a 
aplicação do algoritmo F.F.T.. 


2.2.3 


É altura de passar às transformações expressas 
em (14). Em primeiro lugar é necessário deter- 


minar KK o que em Fortran é realizado através 
das instruções 


XN ==N 
KK = LOG (XN)/LOG (2.) + 0.01 


A forma como foi delineado o programa leva- 
-nos a introduzir uma função indicatriz para 
diferenciar se se está a calcular o conjunto A a 
partir de X ou o conjunto X a partir de A. Se em 
Fortran designarmos a função indicatriz por IND 
então assentamos que IND = 1 indica que estamos 
a calcular X a partir de A e que IND =2 indica 
que se está a calcular A a partir de X. Como já 
foi dito, o primeiro cálculo é feito de A para X, 
logo deve estabelecer-se inicialmente 


IND = 1 
Como em (14) se indica que j é 11,2,3,4 4 KK] 
é necessário recorrer a uma instrução Fortran que 
itere os cálculos. Sendo )J- — j então vem 


DO 70 ] = 1, KK 


O valor 2! entra em (14); vamos designá-lo por 
K2 em Fortran, isto é, 


K2=2.** )] + 0.01 


Vamos admitir que se calculam as partes reais 
em primeiro lugar e depois se repete tudo para 
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"Puração de barrilete com cabeça de furar tripla, com variação contínua de 
entre eixos, para uma das Caldeiras da Central da Refinaria do Norte (Sacor) 


PONTES ROLANTES, GUINDASTES E 
APAR. DE ELEVAÇÃO ESPECIAIS Projecto e fabrico 


TURBINAS HIDRÁULICAS — -——— Fabrico segundo licença de A.C. M. de Vevey, S. 4. 
TURBINAS A VAPOR ——— — ——— Fabrico segundo licença de Brown Boveri, Cie. 
CALDEIRAS A VAPOR ——— Projecto e fabrico segundo licença de Foster 


Wheeler, Co. 
EQUIPAMENTOS E INSTALAÇÕES 


INDUSTRIAIS 


CONSTRUÇÕES METALOMECANICAS MAGUE s..-. 


ALVERCA DO RIBATEJO- PORTUGAL 


TECNICA XII 


Bs q WA 
lj m a ; pi mm nim) El Im 
Ig ai pe 


4 bt Sig pá 


| | o att, Pat It N RR naçio E gnt qnt aim" un mto 
es a ar! mm : À MOI ] Tm Tm l Con MM 1 o] 


tá E TA 


tm 


o! | 
| hit! | I) RE | fd Abba Im lu Mi | 


mt) Nin Eh, mp I 
| ci dl MI 
Fr E: a “ nd Am mm 


É, 


H a apt? | 
Hi it H pr mm , | 
im PN HI 1 7 N' q 
d Fa P IP RRIE 4 Fl 
mms, Jr e, “in 
11 E Fi n !| H ' , 
! Haiti lh , ) À 
a : nb, hi E 


iu 


Gil 
Par Mirim No 
E | “a | 7 mg fi ) & E - 


, ZE a a | rd 


pesei 
] aa 
sá ai Un (| ”, meio “a ii | 
se oii duto O ge O De a 


co cesastdiacemtamas me 


SEPSA - Instalações fabris, Leça do Baílio. 
Vista parcial da Nave de Caldeiraria Pesada. 


ALTERNADORES (de potência superior a 2000 kVA), COMPENSADORES 
SÍNCRONOS, MOTORES DE GRANDE POTÊNCIA (fabricação parcial). 


EQUIPAMENTOS PARA INDÚSTRIAS QUÍMICAS E PETROQUÍMICAS 


Tubagens, tremonhas, reactores, permutadores, cisternas, etc. 


EQUIPAMENTOS PARA INDÚSTRIAS PETROLÍFERAS 


Torres, balões, permutadores, filtros, aerorefrigerantes, tubagens, etc. 


> 
> 
> 
Db EQUIPAMENTOS PARA INDUSTRIAS CIMENTEIRAS 
é 
> 


Fornos, moinhos, separadores, arrefecedores, etc. 


RESERVATÓRIOS 


De tecto fixo e tecto flutuante, de alta e baixa pressão. Para gases liquefeitos. Gasómetros 


EQUIPAMENTOS DE MANUTENÇÃO E TRANSPORTE 


Elevadores de alcatruzes, elevadores de palhetus, transportadores de tela metálica « de 
borracha, contentores, etc. 


ESTRUTURAS METALICAS 


Pontes, cimbres, coberturas, pisos, edifícios industriais, etc. 


MONTAGENS PESADAS 


SEPS SOCIEDADE DE CONSTRUÇÕES ELECTRO-MECÂNICAS, S.A.R.L. 


APARTADO 8 + LEÇA do BAILIO+-TELEX SEPSA P 2616 


TECNICA XIV 


o cálculo das partes imaginárias, ô que em For- 
tran pode ser indicado por 


DO 80KA=1,2 


sendo KA o índice que diferencia na represen- 
tação de números complexos a parte real da parte 
imaginária. 

Outro valor que intervém nas expressões (14) 


Es 


é o de que vamos designar em Fortran por M1 


M1 =N/K2 


de tal forma que as variações dos índices k e L 
em (14) se podem escrever como 


DO 80 [=1, M1 
DO 80 L=1,K2 


se for I--k e L-->L. Se J1I--r então a 
expressão r=(k — 1) 2' + L apresenta-se como 


Ji=kK2*t1I=-K2+L 


Ainda de (14) se conclui que u e v só se defi- 
nem depois de verificar seL<2'7 gu Ls, 
Sendo u ——»] 2 v—-]3 obtém-se 


IFIL—K 2/2) 91, 91, 92 


91]J2=]1 
J3=]J1+K2/2 
GOTO 93 

922]2=]J1—-K 2/2 
J3=]1 

DD us sesmaii 


Em (14) exige-se ainda o conhecimento de 

W(s 281). No cálculo do argumento sabemos 
KK — j N Coni 

que s=L-—l1 e que 2 = “a) * “onjunta- 


mente com (15) é possível escrever 


93 BI=(L - 1)*(N/K2) 
B1 = — 2.*3.1415926* BI/XN 
W(1) = COS (B1) 
W(2 = SIN (B1) 


O problema que temos a resolver agora é o da 
expressão principal em (14) que sendo complexa 
tem que se desdobrar. A forma da expressão é 


(d ,d)=(c,c)+(a,,a)).(b,,b) (25) 
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e, portanto, por separação da parte real e da 
parte imaginária 


d=c+a,.b—a,.b, 
d=c,+a.Db,+a.D (26) 
Verifica-se que o cálculo da parte real e da parte 
imaginária diferem no que diz respeito ao sinal 
da última parcela e à posição de b, e b,. 

Se quando for KA==1 se tiver KI=2 e 
a= +a, e quando KA==2 se tiver KI =1 
ea=-—a,, então as expressões de (26) resu- 
mem-se a 


deaa=Cata,.bea—a.by (27) 


o que corresponde, em Fortran, a fazer 


GO TO(81,82),KA 


81 K1l ==2 
XX = X()3,2) 
AA = A (J3,2) 
GO TO 88 

82 KL=1 
AA = — A(J3,2) 
XX = — X(J3,2) 


88 GO TO (94,95), IND 
94 X(JLKA) = A(JZ,KA)+ A()J3,1)*W KA)— 


— AA*WIK1) 

GO TO 80 
95 A(JLKA)=X(JZ,KA) + X(]3,1)*WIKA)— 
— XX *WIiK1) 


80 CONTINUE 


em que não se esqueceu de indicar pela instrução 
GO TO (94,95), IND se se trata de calcular A 
a partir de X ou vice-versa, 

Após o cálculo de todos os valores de A ou de 
X é necessário trocar o valor da função indi- 
catriz o que se consegue por 


GO TO (96,97), IND 


96 IND=2 
GO TO 70 
97 IND=1 


70 CONTINUE 


Como as instruções de iteração fazem com que 
todos os valores de interesse de ]J, | Le KA 
sejam considerados o programa nesta parte 
está terminado. 


2.2.4 


Vamos nesta fase tratar da saída dos resultados. 
Começaremos pela saída da captação da tabela; 
na instrução Fortran, como não há variáveis a 
indicar, basta escrever 


PRINT 400 


A tabela propriamente dita deve ser dos valores 
de A ou dos valores de X, conforme KK seja par 
ou ímpar, o que se pode determinar pelo valor 
da função indicatriz IND. A tabulação dos valores 
corresponde a uma iteração e, portanto, 


GO TO (98,99), IND 
98 DO 110 ]J1=1,N 
110 PRINT 300,]1,A ()]1,1), A (J1,2) 
GO TO 999 
99 DO 120 ]J1-1,N 
120 PRINT 300, )1,X ()J1,1), X ()J1,2) 
GO TO 999 


com que termina o programa. Verifica-se que 
pelas instruções GO TO 999 se teve em mente 
a possibilidade de utilizar este programa repeti- 
damente: após a saída dos resultados o com- 


putador pode receber novos dados. 
Esta é a última fase do programa o qual se 


encontra completo na Fig. 3. 


2.3 Adaptação para N = 256 


O problema que se enfrenta é a falta de capa- 
cidade de memória para executar o programa 


- 6600 DIMENSION A(128,2),x(128,2),W(2) 
- 6600 999 READ 11,N 
-b624 11 FORMAT(I4) 
-b646 D612 I=1,N 
- 6659 12 READ 13,X(1,1),X(1,2) 
“6778 13 FORMAT(6X,E1L,8,1X,E14,8) 
-BB20 C e q DE DADOS 
-6I20 
6332 A + 
-hILl A a 
b456 DO 2 JeZ.N 
«6363 KKe1 
- B390 | Kiel, **KK+,01 
- 6940 MisN/K2+1 
-6989 DO 3 KAs1,M] 
- 7000 Le J-KI*KA+HKZ/2-1 
- 7104 IF (LL) 5,4,3 
- 7164 3 CONTINUE 
- 7200 5 KKeKK +41 
-7236 e To 1 
-7244 4 Imi-(K2-3)*(N/K2) 
-7352 A(J,2)=X + od) 
- 71424 2 A(J.1 =X(1,1) 
-71532 €C CALCULO Das TRANSFORMAÇÕES 
-7532 KXNeli 
- 7566 KKeLOG(XN)/L06(2,)+0,01 
- 7652 END=1 
- 7676 K2e2,**J+0,01 
-7736 DO BO KA=1,2 
-7748 MI=N/K2 
- 7784 DO HO l=1,M1 
-7796 DO BO Lel,K2 
- 7808 J1=KZTI-KZ+AL 
- 7868 IF (L-K2/2) 91,91,92 
«7960 91 Jladi 
- 7972 JIeJ1+K2/2 


Fig. 3 


10 


anterior para o valor de N imediatamente a 
seguir: 2º = 256. 

Nesta adaptação a principal alteração em rela- 
ção ao programa da Fig. 3 foi a divisão do pro- 
grama em dois. No primeiro é feita a entrada dos 
dados e a sua ordenação. O segundo programa 
refere-se às transformações expressas em (14) e 
à saída de resultados. É introduzido na máquina 
sem se limpar a memória de modo a que os 
resultados do primeiro programa sejam aprovei- 
tados como dados para a execução do segundo. 
Isto consegue-se dando uma instrução DIMEN- 
SION precisamente igual à do primeiro programa; 
assim as posições de memória ocupadas pelos 
resultados do primeiro programa vão ser as posi- 
ções de memória designadas para os dados do 
segundo programa. 

Outra alteração consistiu em substituir a 
variável N pela constante 256 e assim o número 
KK passa a ser constante e igual a 8, o que per- 
mite excluir a subrotina LOG ( ) que ocupa 568 
dígitos. 

Ainda para poupar capacidade de memória, 
embora aumentando o tempo de cálculo, passou 
a haver apenas uma expressão que calcula X a 
partir de A, sendo em seguida passados de X 
para A todos os valores, de modo a no passo 
seguinte ser novamente calculado X a partir de 
A. Como consequência, passou a haver um único 
ciclo de saída de resultados e deixou de ser neces- 
sária a variável IND em todo o programa. 


-ROZO ae 70 33 
- 3028 92 JZweJi-K2/2 
- BO 88 J3w J1 
- 3100 93 Bis(L-1)*(N/K2) 
-8196 Bil=-2,.*3,1415926*81/XN 
-R263 A da o 
8292 2)eSIN(B1 
-3316 60 TO (91,32) KA 
- 3392 AY Kiwz 
— 404 XX X(J3,4 
— 3440 Ah= A(J3,2 
- Bl 76 o 70 3: 
— Bl 82 Klel 
— 8496 Abs-h(J3,2 
-8556 XXe-X(J3,2 
-B616 38 G6 TO (94,95 
E 94 ASR cAjnátia: A Daht IS: VW)*W(KA)-ARMW(K1) 
-899 
- 9000 95 A(JI oe ao KAJ+X( JI,I)MW(KA)=XXAW(KI) 
-9300 Bo CONTIN 
-9405 ne Te te: 97), IND 
— Gu Bh 96 INDeZ 
- 9496 qe 70 
- 9504 97 INDe1 
-9516 70 CONTINUE 
-9552 PRINT 400 
-9564 400 FORMAT RREO PA ER INN DD ROTINA NAO 
-9676 Ge TC (95, 294 1 
-9752 98 DO 110 Jl=l, 
-976h4 HO PRINT 300, J1 EVITA 1),A(J1,2) 
-9296 Go TO 3939 
- 9904 99 DO 120 JI=l, 
- 9915 120 PRINT 300, J1 TATO 1),X(J1,2) 
J0048 ço 78 999 
J00565 300 FORMAT(SX,15,10X,E12,6,10X,E12.6) 


ENO 
PROCESSING COMPLETE 
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Ainda para economizar mais algumas posições 
de memória foram abolidas as variáveis AA e 
XX sendo o papel por elas representado entregue 
a B1 que intervém no programa com duas atri- 
buições distintas, mas que não se misturam. 

Em conclusão, suprimiram-se algumas variáveis: 
como N, IND, AA, XX e algumas instruções 
como, por exemplo, READ 11,N PRINT 400 
e os respectivos formatos. 

A supressão do índice de ordem do argumento 
na impressão da tabela foi feito com o fim de 
reduzir o tempo de impressão da tabela que, 
devido ao elevado número de pontos, começa a 
ser problema. 

Dado que este programa, para além das altera- 
ções já indicadas, é formalmente idêntico ao da 
Fig. 3 só nos resta apresentá-lo: a primeira parte 
encontra-se na Fig. 4 e a segunda na Fig. 5. 


=6600 Ç GRDENAÇAS D 
00 rá (2960 2) .2(296, 2),W(2) 
-6600 
-6612 K 1? 1 pi 
o 12 X 2 A 
s2 
edor fi adia 11 
-6816 2) =X 23 
-6828 
-68h0 =) : 
«68852 | KZe2,**KK4,01 
«6912 His2 S6(K2+ 
-6960 DO 3 M1 
-6972 LosoRINKASKZ/2- 1 
- 7080 AR L) 5,4 
-7136 3 con 
-7172 b) rr 
-7208 6 701 
-7216 4 Imi-(KZ- 1*(256/k2) 
-732h A(J,2)=X a 
- 7396 2 A(J,1)=X(1,1 
- 7504 $ 
=1516 
Fig. 4 
- 6600 prNEuTAON A(256,2),X(256,2),W(2) 
-s600 N=256 
-6636 DO JO J=1,8 
-66h8 K2e2,** 340,0] 
-6708 DO BO KA=1,2 
«6720 MI=256/K2 
-6756 DO 80 I=1,M1 
-6768 DO 80 L=1,k2 
=-6780 J1sKkZ%I=K24L 
-6L4O E RinR2 72) 91,91,92 
-6932 91 
69h Sind aa? 
-6992 OTTO 9 
-7000 92 JZwJ1-K2/2 
7060 J3e J] 
“7072 93 Bi=(L-1)*(256/K2) 
-716t Bl=-2,*3,1415926*B1/XN 
«7240 W Digpio B1) 
-726h W(Z)=SIN(B1) 
-728L Go TO (81,82),KA 
-7364 e! Bi=A(J3,2)*W(2) 
-Th4Zh TO 30 
-7432 8z2 Bi=-A(J3,2)*W(1) 
-7504 80 x(J1 NASA, KAJ+A(J3,I)MW(KA)-BI 
-7840 DO 70 |=1,25 
-7852 A(I,1)ex(i 
-7924 70 A(1,2)=X(1,2) 
- DO 110 Jli=l 2se 
-B0 BO MO PRINT 13, ata W,A(J1,2) 
-8B200 13 FORMAT(6X,E 4,4, O a Y 
-B242 STOP 
HZ Sh ENO 


Fig. 5 
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3. JESTES 


Os resultados obtidos por utilização dos pro- 
gramas Fortran, acima descritos, vai permitir 
avaliar a excelência do algoritmo sempre que se 
conheça o resultado teórico da transformação de 
Fourier. A comparação dos resultados vai ser 
apresentada, com um certo detalhe, em dois casos: 
o primeiro refere-se à transformada de Fourier de 
um rectângulo e foi realizado pelo programa da 
Fig. 3; o segundo caso refere-se à obtenção da 
transformada inversa de Y (4) cujo resultado é a 
exponéncial multiplicada pela função escalão 
(ou de Heaviside) e em que se recorreu ao pro- 
grama para N = 256. 


3.1 Trasformada directa do rectângulo 


Vamos considerar uma função discreta rectan- 
gular Y (8) com a seguinte definição 


Y(B)=1+j0 para B:(0,1,2,...,30) 
Y(B)=0+jo para [e(31,32,...,127) 


Vamos designar por X (x) = [ Y ([) ) a trans- 
formada discreta de Y (É) que se pode calcular 
teoricamente. Vamos apenas interessar-nos pelo 
módulo da transformada que tem a seguinte 
expressão 


IX()|=31 sen (31 7a) 
317ma 
ou 
31 
sen (* a n) 
|X, | =31 Si (28) 
128 
ne(0,1, 127 ) 
pois «= n À « = ain, 
128 


Os resultados deste cálculo tal como a máquina 
os dá estão na Fig. 6. Para comparar os módulos 
da transformada obtida por F.F.T. com a expres- 
são teórica (28) construiu-se a Tabela I, em que 
se apresentam os primeiros 27 pontos. Ainda no 
sentido de facilitar a comparação se procedeu na 


“ 
—s 


-— 


DSI aU. warn 


ooo] do td ué 
+ Ot m O 


) + E + RR > O E + RR O O 
o» wir mM OD a 


REAL 


-310000E OZ 
.2031936E 072 
.199513E 01 
.513160€ 01 


.930786E O 


«MITÚGIE 4] 
«VIS6IL4E UI 
.204L460E 01 


—  B2III0E OU 


«229904E 01 
«JBSI9ZE 01 


-. B31469E UO 


«I4ITIZE 01 
«177301E 01 


i 
mm 
am 
am 
um 
“ad 
Õ 
m 
Õ 
o 


« 904627E-01 
«109801E 01 


«973629€ 00 
« 104832E-05 


«227570E-01 
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TABELA 1 


Valor teórico 


31,0000 
28,0944 
20,3472 
10,2837 
0,9954 
50132 
6,7171 
4,7588 
0,9936 
2,4220 
39523 
3,2228 
0,9852 
1,4090 
2,7401 
2,4832 
0,9745 
0,8626 
2,0462 
2,0362 
0,9602 
0,5172 
1,5885 
1,7284 
0,9432 
0,2786 
1,2587 


«V24GL TE VO 


-195090E 00 


MAGINARIA 


- SO0000E -50 
«1393629E 02 
.<02574E 02 
«SZ267S1E 01 


- 4323)1E 01 
- «Jococ3E DU 
- s4261E 00 
«.ASZ0DLE 01 
«STo4L4GE 01 
«555570E UU 
.210313E 00 
«21004 1E 0] 
.253715E 01 
«JOTIO7E 00 
.435771E-01 
«134131E 01 
.203902E 01 
«SIN4TOE DU 
«131479E 00 
. 186477E 00 
«ITITIME 01 
.923280E 00 
«125944E 00 
«. 391403E 00 
.13"367E 01 
.980785E 00 
.670912E-01 
«JOST4SE DO 
.107423E 01 
« 100000E 01 
.251073E-01 


Valor F. F. T. 


31,0000 
28,0971 
20,3554 
10,2931 
1,0000 
s,0260 
0,7414 
4,7823 
1,0000 
2 4418 
39922 
3,2623 
1,0000 
14333 
2,7948 
2,5402 
1,0000 
0,8881 
2.1143 
2,1119 
1,0000 
O,5411 
1,6684 
1,8237 
1,0000 
0.2969 
1.3483 


Fig. 7 à representação gráfica da curva teórica 
e dos pontos obtidos experimentalmente através 
do algoritmo de F.F.T.. 

Verificou-se que o tempo de cálculo (com exclu- 
são da entrada e da saída) com o programa da Fig. 3 
é de cerca de 18 minutos, no computador IBM 
1620. 


3.2 Transformada inversa da exponencial 


A função Y(2) obtém-se da transformada 
contínua de Fourier da exponencial 


(a +jz=5)-! 


pelo método das dobragens («folding») de tal 
maneira que resulte um Y (£) circunscrito ao 
intervalo limitado T. Sendo assim, a transformada 
inversa discreta de Y (É), será 


X (a) = 375! v(8) | = À eae 
2 
ou, por ser «=n dÁa= — 


— 256 (29) 


Os valores correspondentes à parte real, a única 
que nos interessa, é apresentada na Tabela II que 
também se encontra transposta para o gráfico da 
Fig. 8 na zona mais significativa da transformada 
inversa (n <. 26). No sentido de se obter um por- 
menor mais elevado, procedeu-se igualmente à 
representação do logaritmo neperiano da curva 
teórica (que é uma recta como se pode concluir 
de (29)) e dos pontos obtidos experimentalmente 
pelo algoritmo de F.F. T. (Fig. 9). 

O cálculo dura, se for efectuado num compu- 
tador IBM 1620, cerca de 15 minutos em relação 
à primeira parte (execução do programa da 
Fig. 4) e cerca de meia hora no que respeita à 
cegunda parte (execução do programa da Fig. 5). 
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TABELA Il 


Valor teórico 


0,5000 
0,3955 
0,3129 
0,2475 
o,1958 
9,1549 
O,1225 
0,0969 
0,0761 
0,0607 
0,0480 
0,0380 
0,0300 
0,0236 
0,0188 
0,0149 
o,0118 
0,0093 
0,0074 
0,0058 
0,0046 
0,0036 
0,0029 
0,0023 
0,0018 
0,0014 
0,001I 


Valor F.F.T. 


0,4994 
0,3953 
0,3126 
0,2475 
9,1957 
0,1547 
0,1226 
0,0967 
0,0768 
0,0604 
0,0481 
0,0379 
0,0300 
0,0238 
o,0186 
O,0150 
o,0116 
0,0094 
0,0072 
0,0058 
0,0046 
0,0035 
0,0030 
0,0021 
0,0020 
0,0013 
0,0012 


Fig. 7 


4. CONCLUSÕES 


Da análise das Tabelas I e II conclui-se que a 
aproximação obtida com o algoritmo de F. F. T. 
é boa pois só quando os valores obtidos são 
pequenos (da ordem de 1º/, do valor na origem) 
é que os erros se tornam importantes. Esta con- 
clusão é compreensível visto o computador come- 
ter erros de arredóndamento durante os cálculos ; 
no computador IBM 1620 as representações numé- 
ricas são feitas com 8 algarismos decimais. 

Para valores até 10º/, do valor na origem 
(valor máximo) o erro relativo obtido é inferior 
a 1º/, e para valores até 1'/, do máximo o erro 
relativo cifra-se em alguns por cento (< 10º/,). 

O programa obtido apresenta-se, portanto, com 
interesse mesmo para aplicações na análise de 
sinais apesar do diminuto número de pontos 
imposto pelo computador IBM 1620. Desta forma 
a intenção posta no trabalho parece amplamente 
satisfeita. 


BIBLIOGRAFIA 


— The fast Fourier transform — E. O. Brigham, R. E. 
Morrow — Spectrum, Vol. 4, N.º 12 pp. 63-70 (Dezem- 
bro 1967); 


— Application of the Fast Fourier Transform to Com- 
putation of Fourier Integrals, Fourier Series and 
Convolution Integrals — J.W.Cooley, P. A. W. Lewis» 
P. D. Welch — JEEE Trans. on Audio and Electroa- 
custics, Vol. AU-15, N.º 2, pp. 79-84 (Junho 1967), 


13 


— (Curva teórica 
e o Pontos experimentais 


— (Curva teórico 


. Pontos experimentais 


Fig. 9 


-— fast Fourier Transforms. For Fun and Profit-W.M. 
Gentleman, G. Sand — Proceedings «Fall Joint Com- 
puter Conference 1966»; 


14 


— What is the Fast Fourier Transform ? — Cochran 
et al. — Proceedings of the IEEE, Vol. 55, n.º 10, 


pp. 1664 - 1674 (Outubro 1967). 


TÉCNICA N.º 391 


FRAPIL 


CONSTRUÇÕES E MONTAGENS ELÉCTRICAS 
SARL 


transformadores de intensidade 


—baixa tensão-— 


núcleo magnético toroidal e perdas extra-reduzidas e vasta gama de utilizações 


caixas de polipropileno antichoque 


Pu 
licença 
SEDE E INSTALAÇÕES FABRIS DELEGAÇÃO NO SUL 
CAIS DE SÃO ROQUE RUA CASTILHO 
APARTADO 20 38, ric.-dt, 
TELEFONE 23071 2 TELEFONE 52528 
AVEIRO LISBOA-1 


042.71:2575 . 


TECNICA XV 


ESTAMPARIAS TINTURARIAS 
LAVANDARIAS 


PARA OS VOSSOS TECIDOS TÉCNICOS: 
TRANSPORTADORES, MANCHONS, LAPPING, WOLFRIES, 
FLANELAS, CANEVAS, ETC,, 


EM FIBRAS NATURAIS OU SINTÉTICAS 
CONSULTEM: 


FANAF EL 


FABRICA NACIONAL DE FELTROS INDUSTRIAIS 
SOCIEDADE LIMITADA 
APART. 9 
TELEFS : 52091 PBX O VAR 
TELEG.: FELTROS PORTUGAL ESTRADA DE S. JOÃO 


ESCRITORIO E FABRICA 


Empresa Witramarina de Sondagens & Fundações, L 


(SONDADORA) 


SONDAGENS geológicas e geotécnicas — CAPTAÇÕES de águas 

subterrâneas — REBAIXAMENTOS do nivel freático — CONSOLIDA- 

ÇÕES E ESTABILIZAÇÃO do solo — INJECÇÕES DE CIMENTO 

e outros produtos — INFRAESIRUTURAS de barragens, fúneis e pontes 
— FUNDAÇÕES de todos os tipos 


Única Empresa da especialidade com sede no Ultramar Português 


DELEGAÇÃO - LISBOA SEDE - LUANDA MOÇAMBIQUE 
R.S.Pedro de Alcântara, 1 Rua Serpa Pinto, 60  L.Marques - C.P.982 - Tel.741121 
Tel.32 4693 e 213256 C.P.2178 - Tel.4226 Beira - C.P. 1002 - Tel. 71497 


<CNICA XVI 


C. D. U. 536.212 


CONDUÇÃO DE CALOR NUMA PLACA 


EM REGIME TRANSIENTE 
(DISTRIBUIÇÃO INICIAL DE TEMPERATURAS NÃO UNIFORME) 


RESUMO 


Apresenta-se um método aproximado que permite deter - 
minar, para o regime transiente, a evolução de temperaturas 
numa placa infinita sujeita no instante inicial a uma 
distribuição de temperaturas arbitrária, É (x),e na fronteira 
por condições aos limites de 3,º espécie (convecção para um 
meio a temperatura constante ), 

O metodo consiste na decomposição do problema noutros 
mais simples cujas soluções têm a particularidade de serem 
independentes de E (x), facto que torne viável o tabelamento 
das soluções parcelares, possibilitando a solução do problema 
para qualquer forma de Ffxz), por simples combinação 
linear dessas mesmus soluções. 

Para realização númerica programaram-se as expressões 


deduzidas em linguagem FORTRAN, 


0 — INTRODUÇÃO 


0.1 — Na grande maioria dos problemas clás- 
sicos de condução de calor em regime transiente 
com distribuições iniciais de temperaturas unifor- 
mes, as expressões representativas das soluções 
apresentam, em geral, uma forma acessível a apli- 
cações numéricas. 

Ao pretender generalizar certas soluções para 
casos com distribuições iniciais de temperaturas 
não uniformes, corre-se o risco de comprometer 
a utilização prática das expressões deduzidas, 
pela complexidade e extensão dos cálculos ineren- 
tes a eventuais execuções numéricas. 

0.2 — O estudo da condução de calor em regime 
transiente numa placa infinita segundo y e z, 
de largura 2L(— L< x< L) e condicionada na 
fronteira por um meio a temperatura constante 
T,, constitui um exemplo concreto do que se 
pretende evidenciar. 
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Eng.” Mecânico 
Assist. do 1. S. T.(*) 


SYNOPSIS 


The author presents a method that permits to determinate 
for the transient state, the evolution of temperatures in a 
infinite plate, submited in the initial instant to an arbitrary 
distribution vf temperatures F (x) 

And in the frontier, to limit conditions of 3rd, class 
(convection for a medium with a constant temperature), 

The method consists on the decomposition of the problem 
in others which solutions have the property of beeing inde- 
pendents of F (x). That permits to make tables of partiais 
solutions making easy the solution of any form of F (x) 
for a simple lineal combination of the partials solutions, 

For a numerical realisation the expressicns have been 


written in «FORTRAN» 


Assim, para o caso de uma distribuição inicial 
de temperaturas, uniforme, a solução do problema 


apresenta-se acessível a aplicações numéricas, 
como provam, por exemplo, os conhecidos gráfi- 
cos de Heisler [9,10 ], os quais permitem deter- 
minar, com boa aproximação (!) e para qualquer 
instante, a temperatura no interior da placa. 

Por outro lado, para o caso de uma distribuição 
inicial de temperaturas não uniforme, F(x), a 
solução é ainda conhecida [2], apresentando-se 
com a forma 


oo 


0 — qat 
(0.1) vB lA cos(qx)e aa E 


“ |=1 


+ B, sen (5, x) e 


e cat | 


(t) Esta aproximação, boa para tempos grandes, dimi- 
nui rapidamente quando se aproxima do instante t= o. 


(*) Laboratório de Termodinâmica Aplicada — Departamento de Eng.” Mecânica — Instituto Superior 


Técnico. 
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em que 
L 
/ F(x) cos (q, x) d x 
pe 
Ay —— (0.2) 
I cos” (q,x) dx 
E 
e 
L 
[EGO sen (E, x) dx 
= 
B, = E E (0.3) 
[ sen? (=. x) dx 
as: 


No entanto, expressa a solução nesta forma 
geral, qualquer aplicação numérica forçará certa- 
tamente o utilizador a cálculos morosos e com- 
plexos, dado que: 


— a arbitrariedade da função F(x) torna im- 
possível o tabelamento de expressão 0.1, 
excluindo assim a grande vantagem do caso 
anterior, para o qual a simples utilização dos 
gráficos de MHeisler permite uma rápida 
aplicação númérica ; 


— a resolução dos integrais 0.2 e 0.3 (*) pode 
apresentar-se complexa para efeitos de cál- 
culo, complexidade éssa que dependerá, 
Obviamente, da forma de F (x); 


— se pode tornar necessário utilizar um número 
elevado de termos de série 0.1, ao pretender 
reproduzir, com boa aproximação, a distri- 
buição de temperaturas para instantes pró- 
ximos do inicial, agravando-se a morosidade 
dos cálculos com a irregularidade de F (x). 


Foi precisamente com o intuito de contornar 
estas dificuldades que se desenvolveu, no pre- 
sente trabalho, um processo de cálculo capaz de 
reproduzir, com boa aproximação, a solução 0.1, 
proporcionando ao mesmo tempo uma simples 
e rápida utilização nas aplicações numéricas. 


(?) Coeficientes do desenvolvimento em série da solu- 
ção em termos de funções próprias. 
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0.3— O método que se apresenta, baseou-se 
no caracter linear de equação do calor, 


2 
bdÃ == “O ê - (0.4) 
dt dx 
e das condições de fronteira, 
) 6 -— 
K(s )=+h0 para x=+L (0.5) 
x 


as quais, pelo facto de serem homogéneas, permi- 
tem a utilização do método de sobreposição. 

Sendo assim, procurou-se decompor o problema 
numa série de outros mais simples, de modo 
a que, pela simples combinação linear das solu- 
ções parcelares, se obtivesse a resolução do 
problema na sua generalidade e que, por outro 
lado, essas mesmas soluções tivessem a particula- 
ridade de serem independentes da função F(x), 
garantido-se assim a possibilidade do seu tabe- 
lamento. 

No caso presente, o problema geral foi 
decomposto numa série de problemas parcelares, 
idênticos entre si, cada um deles correspondente 
a uma distribuição inicial de temperaturas da 
forma (fig. 1) 


Fig. 1 
9, (x,0) =1, (x) (0.6) 
com 
I, (x) = do LF (x—b) -—r(x—a)] (0.7) 


em que rF(x) representa a função unitária 
de Heaviside [5], “a, e “b, são parâmetros 
ts ÇD|AL=1L<aÇL) € do uma tempe- 
ratura arbitrária. 
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Seja 0, (x,t) a solução que satisfaz simultã- 
neamente a equação do calor 0,4, as condições 
de fronteira 0.5 e a condição 0.6. 

Uma vez conhecida a solução genérica 0, (x,t) 
dos problemas parcelares, em função dos pará- 
metros T, a, b,k,«x,h, Le T,,o problema pode 
ser automaticamente resolvido pelo método de 
sobreposição. 

Assim, a solução do problema pode ser obtida 
aproximadamente (*) pelo somatório 


gixti= VW 


j=1 


(0.8) 


desde que as N soluções 0, (x, t) correspondam 
a uma série de distribuições do tipo 0.6, que, 
no conjunto e depois de sobrepostas, devem 
reproduzir a função F (x), tal como se encontra 
esquematizado na figura 2. 


Fig. 2 


A partir desta forma, o problema fica reduzido, 
portanto, à obtenção da solução º (x,t). 
0.4 — Criado o modelo físico dos problemas 


mms — ——— 


(*) A expressão 0.8 é a solução do problema na medida 
em que satisfaz simultâneamente a equação do calor, 0,4 
as condições de fronteira 0.5 e, aproximadamente, a distri- 
buição inicial de temperaturas F(x); a precisão de cálculo 
será tanto maior quanto menor for a diferença entre os 
parâmetros a e b, ou seja, quanto maior for o número 
de termos do somatório 0.8. 

Em linguagem mais rigorosa, invocando o teorema do 
máximo das equações às derivadas parciais do tipo para- 
bólico, pode provar-se que o erro máximo cometido para 
qualquer t usando O 8 em vez da solução exacta é menor 
ou igual à diferença máxima entre F(x) e Ó(x,t) para 
t=0, o que é, aliás, fisicamente intuitivo. 


TÉCNICA Nº 391 


parcelares, efectivada a respectiva formulação 
matemática, determinadas as soluções correspon- 
dentes e o erro máximo destas mesmas soluções, 
processaram-se, com o auxílio de um computador, 
todos os cálculos necessários para a elaboração 
de uma série de gráficos que, no conjunto, pro- 
porcionam ao leitor um meio cómodo de aplica- 
ção numérica. 


1 — SIMBOLOGIA 


T, T, — temperatura 


T, — temperaturatura do meio ambiente 
To, To — temperatura arbitrária 
O =T—T, 
0 =T,—T, 
do di To Tu 
e um a PE t, — +, 
? jo To Et a 
é (xt) função definida por 0.3 


x, y e z — coordenadas cartezianas 
t — tempo 
— massa específica 


c — calor específico 
K — conductibilidade térmica 
a = LA ( difusibilidade térmica ) 
pc 
h — coeficiente de transmissão de superfície 
L — metade da espressura da placa 
= EE | número de FOURIER) 
h* = e (número de BIOT ) 
x* = 
L 


4, definidos por 0.3 


q, — valores próprios 

x (t*) — função do tempo adimensional 

v (x*) — função própria 

I. (x) — representação genérica das distribuições 
iniciais de temperatura dos problemas 
parcelares. 

F (x) — temperatura na placa para t== o 
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: | definidas por 4.3.2.2 
b 


0,, A;, B,, D,, E, €, — coeficientes 


+ 


ta | definidos por 4.2.2.1.1 
1, 


>  — para designar transformada de Laplace 
c — para designar a inversa da transformada 
[| ] — bibliografia 


( ) — notas à margem 


2 — FORMULAÇÃO DO MODELO 


2.1 — Hipóteses de cálculo 
Admite-se que: 


a) o material da placa é isótropo; 

b) o calor é a única forma de energia permu- 
tada entre as várias regiões da placa ou 
desta com o exterior; 


Nem 


c) as propriedades físicas K, p, e c são cons- 


tantes e independentes da temperatura; 
d) o coeficiente de transmissão de calor h é 


suposto constante e igual nas duas faces 
limites da placa; 


do 


e) a temperatura do meio exterior é suposta 


uniforme e constante no tempo; 


Do 


não existem fontes de calor no interior da 
placa; 


f 


g) a distribuição inicial de temperaturas é 
exclusivamente função de x; 


h 


T— — 


a função F(x) obedece às condições de 
Dirichelet |5]. 


2.2 — Formulação matemática do problema 
parcelar 


2.2.1 — Com base nas premissas anteriores 
pode-se formular matematicamente o problema. 


Seja 4=T—T,. 


Assim, no interior da placa (—L< x<.L) 
e para t > 0, o fenómeno da condução deverá 
satisfazer a equação do calor 


— = a (2.1) 


sujeita às condições de fronteira 


d 6. 
k ( e )=+h9 parax=+tLet>o (2.2) 
dx 


e à condição inicial 0.6. 


O conjunto da equação 2.1 e das condições 
2.2 e O. 6 definem completamente o problema, 
garantindo a unicidade da solução 9 (x,t) [6]. 


2.2.2 Com o intuito de simplificar o problema 
diminuindo o número de variáveis e parâmetros 
em jogo, optou-se pelo uso das seguintes variá- 
veis adimensionais 


” 5. 
gt = — T, = — 

L um 
tf Dimas at m a 

ç? L (23) 
pe — RL -. b 

K L 


Em função das novas variáveis, o sistema atrás 
formulado toma a forma 


o T* 9º T' asi 1 
| | 
—“=— ara (2.4) 
ot ) x** Ê t* 20 
É gd e — + 1 
j meme 
FR ) = + h*T. para E Sô (2.5) 
l | ” == 0 
T, (x*,0)=1, (x*) as (2.6) 
- TEXT CI 
com 
L(x*)=r(x*t>n)—r(x*—m) (2.7) 


3 — RESOLUÇÃO 


3.1 — Para resolução analítica, recorreu-se ao 
método da separação de variáveis. 


Assim, começou-se por admitir que 


T, (x*,t*) = v(t*) 2 (x?) 


(3.1) 


em que 7 é função apenas de variável t* e > uma 
função dependente, exclusivamente, da inicial x*. 
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Derivando 3.1 em ordem a x* e t* e substi- 
tuindo o resultado em 2.4, obtém-se 


= (3.2) 


em que o primeiro membro é só função de x* e 
o segundo de t*,. 

Como x* e tf são variáveis independentes, a 
equação 2.4 só é satisfeita se ambos os membros 
de 3.2 forem iguais a uma constante À, nos 
intervalos 

Ei dinda (3,3) 
>» 0 

Como h* => 0 por definição, demonstra-se 
matematicamente que a constante À deverá ser 
< 0 [2] o que possui um caracter físico evi- 
dente (*). 

Este facto pode exprimir-se fazendo desde já 


j=—85 (3.4) 
em que à, é uma constante arbitrária. 
Da expressão 3.2 resultará pois 
E Ws 
Sã 5, y = 0 (3.5) 
d x* 
e 
d 2 e 
— =—d dt (3.6) 


Integrando a equação 3.6 pode concluir-se 
imediatamente que 


= — 32 +* 
,=€c,e (3.7) 


em que c, é uma constante arbitrária. 


Por outro lado, a equação 3.5 (equação de 
Helmotz) [1] admite como integrais particulares as 
funções (funções próprias [1]) 


v =D, cos (2,x*) (3.8) 


(') As equações 3.10 e 3.11 revelam que se o expoente 
do exponencial fosse positivo, T; (x*, t*) tendia para CO 
quando t* — 09, o que é fisicamente absurdo. 
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ou 


4, ==E, sen (3,x*) (3.9) 
em que D, e E, são constantes arbitrárias. 

Atendendo às expressões 3.1, 3.7, 3.8 e 3.9 
podemos concluir que tanto 


o 


. (3.10) 


T (xtt*f) =A, e; t* cos (0,x*) 


como 


E (*,t)=B, e-55t* sen (9,x*) (3.11) 
são soluções particulares (”) da equação 2.4, sendo 


A,e B, constantes arbitrárias. 


3.2 Adaptação das soluções particulares às 
condições de fronteira 


Para que as soluções particulares T", (x*, t*) 
satisfaçam as condições de fronteira, a constante 
ô, não pode assumir valores arbitrários. 

Assim, substituindo a solução 3.10 nas condi- 
ções 2.5, conclui-se que 9, só pode tomar valores 
que satisfaçam a equação transcendente 


(3.12) 


Seja q, (i=1, 2,3...) a representação gené- 
rica dos valores de 0, que satisfazem 3.12. 

Do mesmo modo, para que a solução 3.11 
verifique as condições 2.5 os valores de 5, 
deverão ser determinados pela equação 


(3.13) 


Seja :, (i=1,2,3....) a representação gené- 
rica dos valores de à, que verificam 3.13. 

Designam-se por valores próprios os valores 
de à, que satisfazem as equações 3.12 e 3.13, 


(*) Como não se impós qualquer condição restritiva ao 
valor de 5, ' A, e B. , podemos afirmar que existe uma infi- 


nidade de soluções do tipo 3.10 e 3.11 que satisfazem a 
equação do calor 2.4; além disso, pelo facto desta última 
equação ser linear homogénea qualquer combinação linear 
destas soluções particulares continua a satisfazer a equação 
do calor. 


As funções 


O E To. 
T; (x*,t*)=A,e “cos (g,x*) (3.14) 
e 

o % 
T; (28,4 )=Be | sen (5, x*) (3.15) 


além de serem soluções particulares da equação 
do calor satisfazem, portanto, as condições de 
fronteira 2.5. 


3.3 Adaptação da solução à condição inicial 


No instante tf ==0, as funções 3.14 e 3,15 
tomam, respectivamente, a forma 


ts 


> (3.16) 


= A cos (g,x*) 


J: 


A (3.17) 


= B sen (:,x*) 
as quais, só por si, não conseguem reproduzir a 
condição 2.6 correspondente a t* ==0. 

No entanto, como o conjunto de formações do 
tipo 3.16 e 3.17 constitui um sistema ortogonal 
completo [4] e como por hipótese a função corres- 
pondente à distribuição inicial satisfaz as con- 
dições de Dirichelet [5], podemos afirmar que 
uma série do tipo 


i== 00 


X A, cos (g,x*) + B, sen (:,x*) 


[= | 


(3.18) 


pode reproduzir [4] a função T; no instante 
inicial. 

Como, por outro lado, as equações 2.4 e 2.5 
são lineares e homogéneas, a solução do problema 


i=CO ,. 


T,(x*,t*) = > | 


| 
| =| 


sen (mq,) — 


- cos (:,) sen (:,) 


“. 
l 


cos (nc) — cos (m :, 


sen (ng) 


) 9 


pode ser representada, por uma série do tipo 


| == 00 - 48 E 
GC t)= 5 [Ae cos (q,x*) + 
j=1 
— e? (+ 
+ Be "sen (=,x*)] (3.19) 


Resta, pois, determinar as constantes A, e B,, 
de tal modo que T; (x,0) satisfaça a condição 2.6. 

Com base nas propriedades das funções 
ortogonais [5] e atendendo ao especificado na 
condição 2.6, podemos calcular as constantes 
A, e B,, respectivamente, pelas expressões 


m 


| cos (q,x*) dx* 


A = m Fanta (3.20) 
| cos? (q,x*) dx* 
e 
I sen (:,x*) dx* 
B=— (3.21) 


| sen? (: x*) dx* 


n 
as quais, depois de integradas, conduzem a 


sen (mqg,) — sen inq,) 


ma DDD 


(3.22) 
cos (q,) sen (q,) + q, 


cos (n:,) — cos (me,) 


(3.23) 
E — cos (:,) sem (:,) 


Entrando com 3.23 e 3.22 em 3.19, obtém-se 
finalmente 


—9* + | 
— cos(qx*) + 


cos (q,) sen (q,) + q, 


(3.24) 


sen (s,x*) | 
| 
, | 
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4 — ANÁLISE DA SOLUÇÃO T* (x*, tt) PARA 
INTERVALOS PRÓXIMOS DE t* = 0 


4.1 Formalmente, nada havia a acrescentar à 
solução representada pela expressão 3.24. No 
entanto, uma análise mais detalhada das suas 
possibilidades como eventual meio de cálculo em 
aplicações numéricas permite-nos concluir que, 
tal como a solução geral 0.1, apresenta o grave 
inconveniente de necessitar de um grande número 
de termos para reproduzir, com boa aproximação, 
a evolução de temperaturas no interior da placa 
(no intervalo de tempo 0< t* < t, com t7 
pequeno). Esta falha é particularmente acentuada, 
no caso presente, pelo facto da função 1, (x*) ser 
descontínua nos pontos x* =n e x! ==m. 

Para contornar a dificuldade, deduziu-se uma 
outra expressão capaz de reproduzir T; (x*, t*), 
no intervalo 0< tf*< t/, utilizando-se para o 
efeito um modelo físico e um método analítico 
diferentes dos utilizados na determinação da 
solução 3.24. 


4.2 — Formulação do novo modelo 


4.2.1 — Para definição deste novo modelo, 
decompôs-se a sua formulação em dois casos 
distintos, correspondendo cada um deles, respec- 
tivamente, a duas formas particulares da função 
Lt). 

Assim, considerou-se agregado ao primeiro 
caso, que por comodidade de exposição designa- 
remos por CN1, uma distribuição inicial 1, (x*), 
comm = +1enF+1 (fig. 3a) e ao segundo 
caso (CN2) uma distribuição do mesmo tipo mas 
com m=1en5t1(oum+tien=-1) 
fig. 3b. 


Fig. 3 


Seja O < t* < t;, um intervalo de tempo 
que precede o instante inicial. 

Se t; for suficientemente pequeno, pode-se 
admitir que: 
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— para o caso CN1 


T' =0 [OCH <t 
o T: para | 1 ii 
3 x* 0 | xt=—1 
— para o caso CN2 
T; 0 E 
/ x | 
fa: aii Eisi 
o x* A 
Po<t<t (4.2) 
x*==— 1 quandom=1 e nft1 
| xt = quandom £+1en=—1 


Sendo assim, para efeitos de determinação de 
evolução da temperatura TY (x*, t*) nos inter- 
valos 0< t*<. to e-il<x*<. 1, torna-se viá- 
vel a substituição da placa infinita por: 


— um corpo infinito (— co < x*< oo) no caso 
CNI; 
— um corpo semi-infinito 
—1<x*<coquandon =-—-1emft1 


ou 


—co<x*<1 quandom=1en£+1 
no caso CN2. 


Esta equivalência de modelos é, obviamente, 
tanto mais correcta quanto menor for o instante 
tie 
Oportunamente, especificaremos o critério esta- 
belecido para a escolha do valor do instante t . 


4.2.2. Definidos os novos modelos físicos, 
passa-se a apresentar as respectivas formulações 
matemáticas. 


4.2.2.1 Para o caso CN1 pretende-se deter- 
minar uma solução equivalente a T; (x*,t*) 
nos intervalos 


o< yx*< co 
— O € x* < 


O <t*<ti 


(4.3) 


Sendo assim, T; (x*,t*) deverá satisfazer a 
equação do calor 


— (4.4) 
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nos intervalos 4.3. Além disso, deverá reproduzir 
a condição inicial 


TS (x*,0)=1,(x*) para — oo < x*< oo (4,5) 
com 
L(t*)=r(x-—nr(x—m) (4.6) 
Como condições auxiliares de integração, 
podem-se usar 
É od sas === O (4.7) 
e 
T;(o,t*) =0 (4.8) 


4.2.2.2 Sejjarft=x*+H1,m=m+1e 
n=n+1 


Para o caso CN2 pretende-se determinar a 
solução T(r*, t*) nos intervalos ("). 


Sendo assim, Ti (r*, tt) deverá satisfazer a 
equação do calor 


dT) 9º TJ 
= d == (4.10) 
d t d [* 
a condição ao limite 
/ d É 
| ) p* )j=h*T; para rf=0 (4.11) 
ar uy 


(*) Por razões de simetria apenas consideramos o caso 
den=—1lem £*1, correspondente ao intervalo 4.9 


e a condição inicial 
Ti(r*,0)=1I,(r*) para 0O<r*< co (4.12) 


com 


I(r)=r(t—n)—r(rf—m,) (4.13) 
Como condição auxiliar de integração, pode-se 
usar 
TS; (P,t*)=0 para r* = 00 (4,14) 
Este último sistema, tal como foi formulado, 
dificulta a respectiva resolução analítica. Para sua 
simplificação, utilizou-se, mais uma vez, o método 
da sobreposição. 
Por comodidade de exposição designaram-se, 
respectivamente, por CN2/1e CN22, os dois 


problemas parcelares em que se decompôs o 
caso CNA. 


4.2.2.1.1 O problema CN2/1 consistiu na 
determinação de uma função que se designou por 
(r*,t*), correspondente à solução represen- 
tativa de distribuição de temperaturas num corpo 
infinito (— o < y*< co) sujeito no instante 
inicial a uma distribuição 


Ta (r*,0)== 1, (r*) para— oco < r*< vo (4,15) 
em que 
[(r*t)=[(7*—n,))—r(g*—m,) (4.16) 


- 2 (fig. 4), 


com n,-m+t1 em, 
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É EA aÃ - mito, E? 


MINES D'ALJUSTREL — Setúbal 


(Cachofarra) — Carreiras de construção de caixotões flutuantes para 
a execução de um cais de carga de minérios, no Sado 
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CONSTRUCOES TECNICAS, LDA. 


PRAÇA DO MUNICÍPIO, 13-3* TELEF. 366506 LISBOA 2 * DELEGAÇÕES: CABO VERDE - ANGOLA + MOÇAMBIQUE 


a e e 
— “a ME 


a Do = ; 


RE 
Es EEE sa to 


SIDERURGIA NACIONAL — Seixal — Estruturas da Instalação de Laminagem (2.º fase) 
— Grandes Naves Fabris. 


as aa E = eo .“w títítíõíõíõíí88 


devendo a função T; satisfazer ainda, a equação 
do calor 


Ta Ta j 
a 4.17 
dt* dr* 


Como condições auxiliares de integração podem- 
-se utilizar 


T; (— 00, t*)=0 (4.18) 


T, (05,t* )==0 (4.19) 
4.2.2.1.2 O problema CN2/2 consistiu na 
determinação de uma outra função que se designou 
por T, (r*,t*), correspondente à solução repre- 
sentativa da distribuição de temperaturas num 
corpo semi-infinito (0< rf< c) sujeito no 
instante inicial a uma distribuição 
T, (1*,0)=1 (4.20) 


devendo a função T; satisfazer ainda a equação 
4,17 nos intervalos 


nissan (4.21) 
Oo <tt<t 
e a condição de fronteira 
9T, 
7 : =h*T, para r*— 0 (4.22) 
r 


Como condição auxiliar de integração, pode-se 
utilizar 

T, (c0,t*)=1 (4.23) 

4.2.2.1.3 A solução do sistema correspondente 


ao caso CN2, poderá ser representada com rigor 
pela função 


T; (1*,t*) = Th (1º, t*) — Tá; (1*,4%) (4.24) 


Se, no intervalo O < tf* <.té ,se verificar 


Tº =0 r* == 0 (4.25) 


E para 
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e nt 


1 para 


A condição 4.26 pode ser satisfeita escolhendo 
para m, um valor suficientemente grande e para 


tY um valor suficientemente pequeno. Este último 


requesito torna-se igualmente necessário para 
a satisfação de 4.25 e 4.27. 


4.3 — Resolução analítica 
4.3.1 — A solução do sistema de equações 


correspondente ao caso CN1, é conhecida [3] 
e apresenta-se com a forma 


z É sado m— x” n— x* 
TT (xº, 1º) = E [et(ia ) eos (5) | 


(4.28) 


em que erf ( )éa conhecida «error function». 


4.3.2 — A solução do sistema correspondente 
ao caso CN2 pode ser determinada, como se viu, 
pela expressão 4.24. 

A resolução analítica fica pois reduzida à 
determinação das funções T. (1*,t*) e Ty (1º, 1%). 


4.3.2.1 — O sistema correspondente ao problema 
CN2/1 é idêntico ao apresentado no parágrafo 
4.2.2.1. 


A solução é portanto análoga a 4.28, bastando 
substituir x* por r* , m por m, e n por n,. 
Por outro dado, como se pretende a solução 
apenas nos intervalos O<.r*< 2e 0< t*<t/ 
A 


torna-se possível fazer erf ( — = ) = 1, desde 
2 +* 


que se escolha um valor de m, suficiente- 


mente grande, como exige, aliás, a condição 4.26. 
Sendo assim, a solução do problema CN21, 
pode ser representada pela expressão 


* 


To (r* EM) asd 1 f E us 
te”, = | — er (ae )| (4.29) 


4.3.2.2 — Para determinação da forma de função 
T$ (r*,t*) usou-se, como ferramenta matemática 
a transformada de Laplace. 

Assim, aplicando esta transformada em relação 
é variável t* na equação 4.17, obtém-se 


+ d' tp 
Stp—T(rto )=— (4.30) 
dr* 
em que 
1 
= (4.31) 
S 
e 
Ea) Sd ES WE) (4.32) 


Entrando com 4.20 em 4.30 e fazendo P'=S, 
obtém-se 


dt, 
3—P'tp=—1 
dr* 


(4.33) 


que é uma equação diferencial ordinária, linear 
de coeficientes constantes. O integral geral da 
equação homogénea conspondente pode repre- 
sentar-se por 


— Pr* P r* 


+C, e (4.34) 


em que €, e €, são constantes. 
Por outro lado, o integral particular de 4,33 
é igual neste caso a 


(4.35) 


A solução de 4.33, pode determinar-se fácil- 
mente a partir de 4.34 e 4,35, 


Assim, será 


(4.36) 


Por sua vez a representação em termos da 
transformada de Laplace de 4.22 e 4.23 é 


dt ps *=0 (4.37) 
dr? =— t, para [ = 5 
e 
4 aa 1 MS a 
tp (00,5) = Cc pasa =s (4.38) 


que permite determinar o valor das constantes 
Cel. 


Substituindo os valores destas constantes na 
expressão 4.36 e simplificando, encontra-se 
facilmente 

na P r* 
1 
p= 2 —h* = (4.39) 
S S(P+ h*) 


Utilizando uma tabela de transformadas de 
Laplace [3], pode obter-se a solução 


a ) -+- exp (h* 1* + 


T; (r*,t*) = erf =— 
p (r*,t*) = er Err 


= t* hº ) X< erfe ( E + h* Ve) (4.40) 


2V t+ 
em que erfc ( )=1--erf( ). 


4,3,2.3 Entrando com 4.40 e 4.29 em 4.24 e 
fazendo r* =x" +1 en,= + 1, pode determi- 
nar-se a solução correspondente ao caso CN2, 
que toma a forma 


x* + 1 


2Vt* 


T; (1,4) == erf ( 


x ente ( LD 4 pe yes ) — 


2Vt* 


——ee me — —— — 
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) +exp | he Ge) et! | x (4.41) 


o-s(85e)] 
2 2Vt* 
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5 SOLUÇÃO NUMÉRICA 


5.1 Nas aplicações numéricas cuja realização 
depende de expressões complexas, o cálculo pode 
tornar-se extremamente moroso, pelo que se con- 
sidera de interesse, nestes casos, a apresentação 
das soluções numa forma acessível a uma utili- 
zação rápida e simples. 

No caso presente, houve a preocupação de 
efectuar tal procedimento em relação às soluções 
obtidas. 

Assim, todos os elementos numéricos obtidos 
a partir das expressões 3.24, 4.28 e 4.41 foram 
centralizadas muma colecção de gráficos (”) que 
no conjunto, constituem um meio cómodo de 
manipulação numérica. A título exemplificativo 
apresenta-se no apendice B um extracto dessa 
colecção. 

A apresentação dos elementos e meios neces- 
sários a este processamento encontra-se esquema- 
tizada nas alíneas seguintes. 


5.2 Para realização numérica, utilizou-se um 
valor de 0.2 para a diferença entre os parâmetros 
m en. Por razões de simetria, apenas foram cal- 
culados os valores de Tt correspondentes aos 
intervalos 


—1.0< n< — 0.2 


(5.1) 
— 0.8< m<. 0.0 


com m=0,2%n. 

5.3 Os cálculos efectuados foram divididos em 
dois grupos: o primeiro, relativo à obtenção dos 
valores de T. a partir das expressões 4.28 e 
4.41; o segundo, referente ao cálculo dos valores 
da mesma grandeza, mas com base na expressão 
3.24. 


7) A colecção completa de gráficos, pode ser solicitada 
ao Laboratório de Termodinâmica Aplicada (Departa- 
mento de Eng." Mecânica do 1. S. T.). 

A solução numérica abrange os seguintes intervalos 


1.000 <, T* 0.001 (com variação contínua) 
0.0001 t” <. 100.0000 (com variação Contínua) 
0.05 hº 90.00 (com 10 valores paramétricos) 
— 1 e E 1 (com 16 valores paramétricos) 
0.8 <m<. 00 (com 5 valores paramétricos) 
1.0 <a (com 5 valores paramétricos) 
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O campo de acção das duas primeiras expres- 
sões, está, como se viu, limitado pelas condições 
4.25, 4.26 e 4.27. 

Estes requesitos, condicionam o emprego 
destas expressões apenas ao intervalo de tempo 
O< t*<t7, com t, relativamente pequeno. 

Havia, pois, que estabelecer um critério para 
escolha do valor do instante t” , critério esse que 
foi definido com base nas seguintes cláusulas. 

Para t* =t', deverá verificar-se: 


a) Tt < 0.001 calculado pela expressão 4.28 


para x*=— 1 
== — 0.6 
n=:— 0,8 


b) T* 


< 0.001 calculado pela expressão 4.29 
para r*=0 
no = 0.2 


c) T,> 0.999 calculado pela expressão 4.29 
para 


T* = 


n, = 0.2 


d) T, > 0.999 calculado pela expressão 4.40 
para 


p* = 2 


Com base neste condicionalismo, optou-se por 
um valor de t, = 0.002. 


5.4 O grande volume de cálculos a realizar 
obrigou a utilizar um computador, sem o auxílio 
do qual era práticamente impossível a solução do 
problema numa forma gráfica. 

Para isso programaram-se, em linguagem «FOR- 
TRAN», as diferentes fases do problema. 


5.5 Para finalizar, apresentam-se, de um modo 
sucinto, as diferentes etapas em que escalonamos 
o processamento numérico. 


5.5.1 Cálculo dos valores de T, , no intervalo 
de tempo 0.0001 < t* < 0,002. 

Utilizaram-se para o efeito as expressões 
4.28 e 4,41. 


5.5.2 — Determinação do número de termos da 
série 3.24, necessários à sua utilização. 


Para resolução desta etapa elaborou-se o pro- 
grama n.º 1 (Quadro 1). As instruções deste 
programa foram concebidas de modo a perguntar 
ao computador qual o número de termos, N, da 
série 3.24, a partir do qual a diferença entre os 
valores de T/ por ela calculados e os correspon- 
dentes à função 1, (x*) passava a ser inferior 
a 0.005 (5). 

A comparação dos resultados foi, evidente- 
mente, efectuada para t* = 0, instante em que 
a distribuição de temperatura é conhecida e 
coincide com 1I,(x*); além disso, utilizou-se 
hf== co, pelo motivo de ser extremamente 
simples, para este valor de x*, o cálculo 
(') de valores próprios (:, e q;) correspondentes- 

A resposta do computador foi N = 55. 


5.5.3 — Cálculo dos valores de :, e q; 

Com base nas equações 3.12 e 3.13 elaborou-se 
o programa n.º 2 (Quadro ll) que permitiu 
a determinação dos primeiros 55 valores de :, q, 
(10), para diferentes valores de h*. 


5.54 — Cálculo dos valores de T; no intervalo 
0.002 <. tf < 800. 


Para realização desta última etapa utilizou-se 
o programa n.º 3 (Quadro III), tendo servido 
de dados os valores de E e q, anteriormente 
cálculados. 


(*') O erro da expressão 3.24, para um valor finito de 
N, é máximo para t” == 0. Imposta uma diferença máxima 
de 0.005 relativa a tf?=0, o erro será inferior para 
t* == 0,002, instante a partir do qual se começa efectiva- 
mente a utilizar a expressão 3.24. 

A comparação realizada entre os valores dados por 3.29 
e 1, (x*) foi efectuada para vários valores de x" & m ou n. 
Para xX*=m ou n (correspondentes a pontos de disconti- 
nuidade da função 1, (xº) ) a série 3.24 apresenta, para 
t* ==o, erros de 50'/,, mesmo para valores grandes de N; 
este facto não compromete de modo algum a precisão de 
cálculo efectuado, dado que, para t” == 0.002 o número 
de termos escolhido pelo computador, segundo o critério 
apresentado, confere a 3.24 uma precisão suficiente em 
comparação com a dos gráficos utilizados na apresentação 
das soluções (Apêndice A). 


(') Para h*f = oo 
g=nr (i=1,23....) 


q = A — t ) = (1=1,2,3.0.) 
É 
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Ao elaborar este programa, houve a preocupa- 
ção de reduzir ao máximo o tempo de operação 
do computador. 

Para isso, as instruções foram elaboradas de 
modo a que fosse utilizado, na série 3.24, apenas 
o número de termos indispensáveis para se manter 
a precisão de cálculo estabelecido para t* = 0.002. 

Para este instante o computador utilizou os 
54 (!) previstos pelo processamento n.º 3. Em 
contrapartida, para tf == 800, utilizou apenas 
um termo, tendo-se verificado, no conjunto dos 
cálculos, uma redução de tempo de 80º/,. 


QUADRO | 
PDG FCRTRAU C? 
START 
- 500 PRINT 1 
-5614 | FURMAT(IZHPRUGRAMA NUMERU 1) 
-6670 DG 500 J=1,19,% 
«bis RJ=J 
-571) X=-1.+.1*RJ 
-6700 Tisl. 
-D77S Y=5.HWiIs9z/2, 
-bu 4 De £20uU 1=),4949 
-b026 Rl=l 
-6922 EsB+Y 
-b95u SQ=SIN(C) 
-6982 Y1= (SO SIN(. q PSA 
-7114 rg dad 
-7162 Es tA 
-7196 CE=COS(E 
-1222 Cl= (cos (. 8*E )-CEJ*SIN(E*X) 
-7342 C2=E-SIN(E)*CE 
- 1402 C3=C1/C2 
- 7438 T1=11+13+€3 
7426 IF(SENSE SWITCH 1)150,130 
- 15006 150 PRINT 2,71 
- 7530 2 FCRMAT(F1O,4) 
- 7552 130 1F(X-.9)15,20,20 
7520 I5 Z1=7T1 
-7632 Z2=hBS(Z1) 
-7555 IF(Z2-,005)10,200,200 
-J 7êh ZO T3=Ti-1. 
7760 Z4=ABS(Z3) 
773, IF(Z4-,005)10,200,200 
-7 551 «LO CONTINUE 
-J5va ju PRINT 14,X,RIL, 
-7936 19 FCRMAT(ZHX= Eh: 1. [ .3HRI=,F4.U,3HT]=,F 7.4) 
-S01U 5uu CONTINUE 
- 8046 STCP 
=4USO END 


(1, Os resultados obtidos podem ser solicitados para 
o Laboratório de Termodinâmica Aplicada (Departamento 
de Eng." Mecânica do 1.5.T.). 

Os valores de :, e q; foram calculados para a seguinte 
gama de h* 
— (1,05 — 01 — 0.2 — 0.3 0.4 — 0,5 
— 09 — 1.0 — 1.2 — 1.4 — 1.6 1.8 — 2.0 — 2,5 — 3.0 — 
— 40 — S0 — 60 — 7.0 — 80 — 9.0 — 10.0 — 120 — 
— 140 — 16.0 — 18.0 — 20.0 — 25.0 — 30.0 — 4,00 — 
— 45.0 — 50.0 — 60,0 — 70.0 — 80.0 — 90.0 — 


0.6 ais 0.7 NEES 0.8 ni 


(!!) Em rigor o processamento n.º 3 previu 55 termos, no 
entanto por razões relativas com a capacidade de memória 
da máquina utilizada, teve de se optar por um número de 
termos ligeiramente inferior, ou seja, 54. 
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QUADRO II QUADRO TI 


PDS FURTRAN C2 START 
START -b6600 L PREGRAIA NUMERO 3 , 
-6600 € PRUGRAMA NUMERO -6600 DIMENSION 0(5h),E(S4 on DRA ra e 
- A5NO DIMENS ICN (55. 5, E(55,5),MI(S) -6600 DIMENSION Y(54),2(54),CQ(4,54),SE(M,5h) 
-660C De 1 J=1,5 =8600 DO 1 L=1 Kg 
«6612 | READ 2.H1(J) -6512 | READ 2,6 
«A696 2 FORMAT(F 4,2) -56 6 2 ERNATIEASS 
5713 19 E 3 Jel, “5713 DC 3 L=1,5 
-5730 Ve/HI(S -6730 8 READ 2 ih 
-5799 PRN 390.4 51h DO 12 K=1,35 
Bulb uy FLRMATI///.2HJ=, 13.1) -A926 12 READ 1,TÍK 
-456L DC 12 1=1,55 -6910 1! FCRAT(F 3.3 
-bUHU Ri=| -6932 DO 13 M=1,5 
5916 IF(SENSE SWITCH 1)120,121 - S9Ll, 13 READ 14 .U(M 
“0930 ZU XI=(RI=-1,)"3, 141590 -7024 14 FORMAT(F4,1) 
3904 IF(SENSE SWITCH 1)142,121 - 7050 DO 15 |=1,4 
« FU UM 12) XI=s(RI-1.)e3. JU IG9Z+3. 4 ISUA/2, - 7062 15 READ 14 x(1) 
«7100 |Zz DE 13 12=1.9999 - 7146 READ 16,H1 
-Jh14 Riz=Iz -717U 16 FCRMAT(FL.2) 
-JIbs X=X1+1 ,“RIZ -7192 DC 17 41,5 
-JIzo IF(SENSE SWITCH 1) Iuul, IuuL -7ZUM4 RN=U(M)-,2 
7/15 1002 Y==(X/H -7264 RH=U(.4) 
= Tt ul S=SIN(X -7300 PUNCH 13,RN,RM 
=Tido C=CUS(X =7336 13 FORMAT(///,6X,3HRN=,F4,1,3X,3HRi=,F4,1,/) 
«7514 TeS/€ -7424 IF(SENSE SWITCH 1) 19,20 
73h!) IF(SENSE SWITCH 1) 1001,125 - TU ll 19 PRINT 1U,RN,RH 
-t30L OU) Y=(X/H) -7480 20 DO 400 L=1,54 
- 740 SeSIH(X) -7492 8Hi=RHSO(L 
-742! C=CLS(X) -7552 QN=RN*Q(L 
-7h452 CCT=C/S =-1612 EM=RHSE(L 
«7h 125 IF(SENSE SWITCH 1)126,127 =7572 ENsRN*E(L 
- 7505 126 1F(Y-CCT)13,16,I15 =7732 AT=SIN(OM)-S IHN(O) 
7576 127 1F(Y-=T)1S,15,13 —7 UNA ed 7? + 
=? 13 CONT IHUE —7943 Y(L)sA1/A2 
«7680 15 XeX- le -8008 8B1=COS(EN)-=COS(E 
7715 50 X=X -3080 SEUS OINEE CA ecos tetu)) 
-775% tF(SEnse SWITCH 1)123,129 =—236 2(L)=B1/82 
-7772 123 Y -R29h 400 CONTINUE 
-7303 S=SIN(X) -3332 pe ap ei A 
-7832 C=Ces(x) 23h44 Du L=1 
- 7856 CCT=C/S -3356 ço r Ki ele | 
- 7902 IF(SENSE SWITCH 1)120.129 2512 ai e pes INR 
-791% 129 Y=-(X/H 266" noz é 
- 7960 S=SINOO 2704 L01 CENT UNIR 
-793h4 C=Custx) = 740 DO 5 Iut,h 
- 3008 1=5/€ Da PUNCH 5 XD), H1 
- 8064 131 1F(Y-CC7)13,16,19 Cao IE(SE SE SU ITEM 1)21,2 
-8134 132 IF to thid. 16,19 3914 Z1 PRINT 6.X(1) H 
-B2U0 13 G6 To = 1873 22 DC 200 Kel, e 
-B208 19 tb -A999 Tl=u. 
-8256 51 X=X+,01 «2002 DL 300 L=1,5h 
-Bey2 IF(SENSE SWITCH 1)133,134 0 Vl EEI=E(L)S*2*T(K) 
-8312 133 Y=(X/H) -9110 IF(EE1-S6.)305,303,309 
-B346 S=SIN(X) -9175 301 EEZ=EXP(EEI) 
-8372 C=C0S(X) -9202 309 IF(EE1-55,)510,310,311 
-8396 CUT=C/S 9270 310 EE=1,/EEZ 
-B432 IF(SENSE SWITCH 1)135,134 -9306 311 IF(EE1-55.)312,312,315 
-B452 134 Ym=(X/H -3374 313 EE=U 
-8500 ii -9336 312 EdI=I(L)SAZATIK) 
-8524 E=CUS(X -9452 IF(E91-55.)301,301,3.2 
-BS LE —9550 301 EC2=EXPIESI) 
-B5 By 135 TetiEnse SWITCH 1)136,137 -9574 302 IE tEA coa. )303,303,3U05 
-96h2 303 E9=1./ES 
-967U 305 ETEÓIAS. vo! 307,305 
- R6O4 136 Is fig pb iA 16,22 =97h5 307 ao bo “EQ 
- 8672 137 IF(Y- -9575 F2=Z(L)*SE(I,L)XEE 
- 3740 21 JOO1Oo TisTI+F1+F2 
-8748 22 Xme,UI+X 
-2796 2 XeX+,003 
| -BH32 IF(SENSE SWITCH 1)138,139 
-9352 139 Y=(X/H) 
«BB SeSIN(X) Juoss 300 CONTINUE 
“ =3912 C=CrS(x) 3009h 306 SeT(!) 
“2936 COTeC/5 J0130 PUNCH 30,S,TI,L 
-4972 Md ste SWITCH 1)141,139 J0172 30 FCRUAT(3X,2HT=,F3,3,3X; qria, F6.4,3X,2HL=,12) 
- 3994 139 Y=-(X JO266 IF(SENSE SWITCH 1)31,20 
-90409 ai ) J0236 31 PRINT 30,S,TI,L 
— SU tl Ce CoS(X J0334 200 CONTINUE 
- 90 35 T=S/C J9370 SG CUNT INUE 
-9124 141 IF(SENSE SWITCH 1) 142,43 JOLO 6 STOP 1000 
-914k 142 sá Y-COT)24, 16,16 MUITA 17 CONTINUE 
=GE lê 14> 1F(Y-=T) To, 16,24 JS STCP 20009 
-92£0 24 GC TO 52 JI466 ED 
-92b0 ló L=Ri 
“9324 IF(SENSE SWITCH 1) 149,150 
= 3h t4g O(L,J)=X 
«9416 IF(SENSE SWITCH 1)151,150 


=9436 150 E(L,J)=X 

-9508 151 IF(SENSE SWITCH 1)152,153 
-452 6 152 PUNCH 29,0(1,9) 

-9612 29 FORMAT(F 3.3) 


9634 PRINT 290,.,J,0(1,9) 
-97h2 290 FORMAT(3X, Ena PR 11H50 12,2H)e,F 6,3) 
-9812 IF(SENSE SWITCH 1)12,153 


-5832 153 PUNCH 30,E(L,4) 

-9916 30 FORMAT(FS.3 

-9936 PRINT 291,L,4,E(L,4) 

30046 291 FORMAT(3X,ZHE(, IZ, TARTE 2H)=,F3.3) 
J0116 12 CONTINUE 


J0152 3 CONTINUE 

JO(RA IF(SENSE SWITCH 1)15h4,155 
J0203 154 STEP 1000 

J0220 GC 760 10 

J0228 155 STOP 2000 

JOZ4O FND 
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ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL **) 


— Breve nota mensal 


As afluências ao conjunto do sistema registadas 
durante o mês de Agosto situaram-se abaixo da média, 
tendo totalizado 97 GWh, valor a que corresponde um 
coeficiente de produtibilidade hidroeléctrica de 0,66 e 
que tem uma probabilidade de ser excedido da ordem 
dos 89º/,. 

Estas afluências abaixo da média resultaram não só 
das reduzidas afluências «internas» mas também do 
facto de o caudal transitado em Miranda se ter situado 
sensivelmente abaixo da média, essencialmente pelo 
reduzido desarmazenamento verificado em Ricobayo, 

As afluências acumuladas desde o início do ano 
hidrológico correspondem a um coeficiente de produ- 
tibilidade hidroeléctrica de 1,35. 


LI — Elementos gerais (GWh) 


a) Mensais Variação 


1968 | 1969 | “o 


Produção hidráulica (Ph) . . 


509,9 | 469,9/+ 52 
Produção térmica (P+)..... | 127,3 12,97 |— 90 
Produção total (PT)... ...... 437,2 | 482,6 + 10 
Energia recebida de empresas 
não pertencentes ao RNC (Er) 0,2 0,8 + 50 
Exportações (Ex) .,.... en dos 1,7 3,0 + 76 
Importações (1)... «cc... ...| 13,6 3,0/— 79 
Saldo importador (Sj)....... + 11,9 0,0 :—100 


Cons. em bomb, hidroelé. Ch) (1) 0,0 0,U O) 
Produção para con- 
sumos perman, 
Produção para con- 
sumos não perman, (Pcnp) ... 
Total «esse gusas CEA qa 
Coeficiente de hidraulicidade 


(Pep) e... 


+ 32 
+ 7 


36,8 | 480 
449,3 | 489,9 
0,78) 0,66 


NOTA: 


(1) O aumento percentual da produção para consumos permanen- 
tes, tendo em conta a incidência dos domingos e dias especiai', 
é respectivamente de 7,6 e 89º. 


Ill — Diagramas de carga dos dias característicos 


4.º feira: 


14683 1969 
Produção hidráulica (DP )MWb [12879 [16781 
Produção térmica (Po)MWb | 3050 TO 
Produção total (Pr)mwhb [15409 16 851 
Trocas com Í Export.  (Ex)MWh 30 65 
Espanha ÀImport. ( 1)MWh | 380 | 126 
Cons. em bomb, hidroel. (Ch) MWh 4 U 
Prod. p.* cons, perm.  (Pep) MWb| 14951 15 583 
Prod.p.!cons.não perm.(Penp)MW»s| 1228 1348 
TOTAL . Pr+-(-Ex) MWh [197.59 16 411 
Potênciamáx. MW| “42 911 
“e Potência min. MW 434 462 
I- 
Ê : Etr Utiliz. da ponta horas 18,7 18,6 
z E Factor de carga 0,78 0,77 
ta Potência máx. MW, 793 864 
E p Potência min, MW 3º0 404 
eo SP Utiliz. da ponta horas 18,3 18,0 
Factor de carga 0,76 0,75 


AGOSTO 


Albufeiras: 


Alto Rabagão +... ... 
CANADA é cimo = o é é 
VÓRUA NONE » no vo ms 
Salamonde . . . 0. q 
Canada cc snecscs 
MED mu areas q té DAS! FTA IR É 
CRB. eco ua css à a 
Castelo do Bode +... ... 
CM = dm Ea a 


Lagoa Comprida +. +... «|(M3B18 

Manta LUZIA « «mis «és 51,0 82,8 
Pracana » co co o. er. 9,4 72,9 
DOM a ed En ms Ee (3) 9,0 62,5 
E com A, Rabagão . . .| 1753,9 52,6 
Fotal sem A, Rabagão. . 903,8 18,6 


NOTAS 


(1) Coeficiente de enchimento em energia das albufeiras, 


(2) Inclui 3,2 GWh armazenados em Vale do Rossim no início do 
mês e 2,3 GWh no fim do mês, 


(3) Inclui 2,5 GWh armazenados no açude do Poio no início do mês 
e 2,2 GWh no fim do mês. 


(*) Elementos extraidos das estatisticas mensais do Repar- 
tidor Nacional de Cargas (R. N, C.). As produções e 
os consumos das empresas do R,N.C. representam 
cerca de 949/, dos totais do País, 
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DETERMINAÇÃO DAS CARACTERÍSTICAS DE CORTE 


DE ROCHAS É DE DESLIZAMENTO DE DIACLASES “ 
INFLUÊNCIA DA OBLIQUIDADE DA SUPERFÍCIE DE ENSAIO 


RESUMO 


Descreve-se a aparelhagem e o método de ensaio uti- 
lisado pelo LNEC na determinação das caracteristicas 
de corte de rocha e de deslizamento de diaclases em 
provetes prismáticos com a dimensão aproximada de 
14X 18x 18 cmº. 

Estuda-se teôricamente a influência da obliquidade 
da superfície de rotura ou de deslizamento nos valores 
da coesão e do àângulo de atrito e analisam-se alguns 
casos particulares, 

Verifica-se, por meio de ensaios, a boa concordância 
existente entre a teoria e os resultados experimentais. 


1—- CONSIDERAÇÕES GERAIS 


O Laboratório Nacional de Engenharia Civil 
(LNEC) possui no Serviço de Barragens (SB), 
Divisão de Fundações e Túneis (DFT), uma 
máquina de ensaios que permite a determinação 
das características de corte de rochas ou de des- 
lizamento de diaclases a partir de ensaios reali- 
zados em Pprovetes prismáticos com a altura 
máxima de 14 cm e as dimensões máximas da 
secção transversal de 18 x 18 cmº, 

A máquina de ensaios (fig. 1) é constituída 
por um quadro metálico (A) com cerca de Im 
de lado. Aplicado à travessa superior existe um 
macaco hidráulico (B) com a capacidade de 40 
toneladas ligado a um sistema de manutenção 
da carga constante (C) e aplicado a um dos 
montantes um outro macaco hidráulico (D) com 
a capacidade de 50 toneladas ; ambos os macacos 
hidráulicos transmitem as forças ao provete de 
ensaio por meio de rótulas (E) de forma a permi- 
tirem uma distribuição quanto possível uniforme 
da tensão. Uma bomba hidráulica (F) ligada aos 
dois macacos permite, por meio de um comuta- 
dor (G), aplicar forças a cada um deles, separa- 


por  F. PERES RODRIGUES * 


SYNOPSIS 


The testing equipment and methods used by the 
LNEC for determining the shear throughout the rock 
and the slide along joints in prismatic specimens about 
14X 18x 18 cm” are described. 

The influence of the obliquity of the failure or sliding 
surfaces on the values of the cohesion and the angles of 
friction is theoreticálly studied and some particular 
cases are analyzed, 

According to the tests carried out, the agreement 
between theoretical and experimental results is good. 


damente. O provete de ensaio é colocado no 
interior de dois pequenos quadros metálicos (H) 
rígidos, o inferior fixo e o superior móvel, 


Fig 1 


deixando entre eles um afastamento julgado con- 
veniente para o ensaio, em geral 1 a 3cm, por 


(!) Trabalho apresentado em inglês ao Simpósio Internacional de Mecânica das Rochas, Madrid, Outubro de 1968 
e em português às 2.º Jornadas de Engenharia e Arquitectura do Ultramar, Luanda, Maio de 1969, 
* Engenheiro Especialista, Chefe da Divisão de Fundações e Túneis do Laboratório Nacional de Engenharia Civil, 


Lisboa. 
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